*  V1  - 

I  ...  4 

j 


MUh  U  .  JU 


Heavd-uced 
tuf-  ttte 


!  f  i:  •■ids  technical  information  agency 

ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 


ts.i 


PTCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  my  purpose 
other  than  in  connection  vith  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  ary 
obligation  whatsoever;  ana  the  fact  that  the  Govern¬ 
ment  may  have  foredated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


WADD  TECHNICAL  REPORT  (1-133 
VOLUME  II 


C'3-  -  /- 


.ox 

tj 


PASSIVE  AERODYNAMIC  ATTITUDE  STABILIZATION 
OF  NEAR  EARTH  SATELLITES 

Volume  II 

AERODYNAMIC  ANALYSIS 


PAUL  H.  DAVISON 

NORTH  AMERICAN  AVIATION,  INC. 
COLUMBUS ,  OHIO 


JULY  1961 


uzii 


AERONAUTICAL  SYSTEMS  DIVISION 


Best  Available  Copy 


NOTICES 


Vi  her  Go 'eminent  drawings,  sp*‘cific.iticr  s,  or  oihei  asta  are  used  tor  ^ny  purpose 
ether  than  i"  correction  with  i  del'. many  rela’ed  Government  procurement  operation,  the 
United  States  CJo/ernment  thereby  incurs  no  responsibility  nor  any  obligation  whatsoever, 
and  *m*  fact  that  the  Government  may  have  formulated,  furnished,  c-r  ;n  *ny  way  supplied 
me  said  dra«virgs,  specifications,  or  other  data,  is  not  to  be  regarded  by  implication  cr 
otherwise  as  in  any  mam^er  licensing  tr  e  holder  or  any  other  person  or  corporation,  or 
conveying  any  »*ig  its  or  permission  to  manufacture,  use,  or  sell  any  patented  invention 
that  may  in  any  way  be  related  thereto. 


Qualified  requesters  may  obtain  copies  of  this  report  from  the  Armed  Services  Tech¬ 
nical  Information  Age  icy,  (ASTI A),  Arlington  Hall  Station,  Arlington  12,  Virginia. 


This  report  has  been  released  to  the  Office  of  Technical  Services,  U.  S.  Department 
of  Commerce  WaJtoi  gton  25,  L>.  C„  for  sale  to  the  general  public. 


Copies  oj  ASD  Technical  Reports  an  i  Technical  Notes  should  not  be  returned  to  the 
Aerc nautical  Systems  Division  unless  return  is  required  by  secur.iy  considerations,  con¬ 
tractual  obligations,  or  notice  on  a  specific  documen*. 


\V\DD  TECHNICAL  REPORT  61-133 
VOLUME  II 


PASSIVE  AERODYNAMIC  ATTITUDE  STABILIZATION 
OF  NEAR  EARTH  SATELLITES 

Volume  11 

AERODYNAMIC  ANALYSIS 


PAU L  11  DWJSON 

MW Til  AMERICAS'  AVIATION,  JSC 
COLUMBUS,  OHIO 


JULY  mi 


FLIGHT  DYNAMICS  LABORATORY 
CONTRACT  X.  AF  33(616)-71U9 
PROJECT  Nr  1366 
TASK  Nr  13967 


AERONAUTICAL  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIK  FORCE 
'VKICHT  PATTERSON  AIP.  FORCE  BASE,  OHIO 


V^Cregor  0  Werner,  Lx:. ,  Ga/tan,  O. 
■*00  -  November  135X  -  ’u-341 


Tne  research  reported  herein  was  cerioxr^c  ty  ;.or*r 
/merjran  Aviation,  Inc. ,  Columbus ,  Chic  for  tne  Hyp^.-ic'-'c 
rli^ft  Se»’t  rlirr.t  Branch.  of  the  Fligst  Dynamics  Labor  ,1- 
tor;/,  Wrirht  Air  Devclorr.ent  Oit'ision.  Tne  work  was  accom¬ 
plished  unaer  Air  Force  Cortract  No.  Ar  33(6i6)-71CO, 

Project  No.  1  '66.  Task  No.  13967,  4A  Study  of  Aorodynar.icaily 
Oriented  and  Staoilized  Satellites."  This  research  va> 
carrieo  out  cy  the  Engineering  Research  and  Aemtr  errodyraTjcs 
Development  Groups  of  the  Columbus  Division,  Kor‘h  kr>'w'ir?.r 
Aviation,  Inc.,  with  Dr.  D.  K.  Schrcllo,  Engineer  irg  nesearch 
Grouc,  is  Project  Engineer.  Mr.  Joseph  Ordrejka,  flight 
Dymurjcs  Laboratory,  was  WADI)  Project  Engineer. 

The  results  of  this  stony  are  reported  in  a  series  ff 
three  voluiues  cf  which  this  is  Yolune  IJ.,  The  other  rerc’ts 
in  this  series  are; 

VOLUME  I:  "Lihrationc  due  to  Combined  Acr\/d/narsic 

and  Gravitational  Torques,"  oy  D.  K. 

Scmellc 

VOLUME  III;  "Matnematical  Techniques  and  Ccspjter 
Program,"  by  0.  C.  Juelich 

Particular  credit  snceld  be  given  to  Mr.  J.  E.  Shuter 
for  his  contributions  to  tne  analysis  of  the  more  complex 
satellite  sr.apes  presen* ed  in  this  report. 


V4DD  T'  61-133  Yol II 


ii 


ABSTRACT 


Acrodi/r.amic  drag  ai  i  stability  coefficients  ol  several 
simple  s»atelli*«>  configurations  are  calculated  using  an  approxi¬ 
mate,  but  nearly  exact.'  iree  molecule  OiPCi'j  .  The  rootle  1  < 
examined  are  che  cylinder  sphere,  pointed  cone,  truncated  cone, 
power  law  fcoay hemisphere* cylinder  and  r one-cylinder,  From  the 
viewpoint  of  .laxtmn  aerodynamic  effects.  favoi<ible  center-of- 
mss  locations  are  oo tamed  by  concentrating  a  Lcrugeneous 
pay lead  forward  of  the  center-of-volune.  For  calcula*ion  pur¬ 
poses,  the  aft  portion  of  the  models  are  assumed  to  be  rigid, 
but  massless,  shells  cr  skirts. 

Basest  upon  the  analysis  of  satellite  angulat  motion,  the 
aeroaynasie  pitching  moment  slot>e  is  combined  with  the  pitching 
r.orent  of  inertia  to  form  a  stability  parameter  which,  when 
uxiniced,  acscnt.es  the  optimum  external  geometry  of  each  model 
for  a  specified  payload  mass  distribution.  The  various  optimum 
configurations  are  shown  to  have  low  fineness  ratios  (order  of 
unity),  therecy  permitting  the  gravitational  torque  to  augment 
the  aerodynamic  torque.  Among  the  configur itions  analyzed,  the 
pointed  and  truncated  conical  bodies  and  th?  power  lav  bodies 
are  the  most  stable. 

The  aerodynamic  d&rping- in-pitch  derivative  is  also  calcu¬ 
latea  for  the  cylinder,  sphere  and  cone.  Aerodynamic  damping, 
however,  has  an  insignificant  contribution  to  the  satellite's 
angular  motion. 
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INTRODUCTION 


Although  Me  effects  of  rarefied  aerodyna-dc  forces  on  satellite 
p*rf onaar <*e  (c .sr- ,  orbit  decay)  have,  by  necessity,  been  extensively 
considered,  there  has  been  litt.e  concern  over  the  utilisation  of  these 
forces  fer  passive  attitude  control  and  stability  of  satellite  vehicles. 
This  rather  raid  interest  appears  to  have  been  primarily  due  to  the 
uncertainty  of  the  characteristics  of  the  upper  atmosphere  (say  in  excess 
of  100  miles  altitude)  as  veil  as  to  the  reluctance  to  accept  relatively 
unproven  (or  unknown)  methods  of  estimating  rarefied  aerodynamic  charac¬ 
teristics.  Kith  improved  and  more  frequent  experimental  investigations, 
however,  the  physics  of  the  upper  atmosphere  is  becoming  more  clearly 
understood.  Likewise,  tnreugh  the  inspiration  of  Schaaf,  Patterson, 
Stalder  and  their  colleagues  (Refs.  1,  3,  5,  9,  9,  12,  13.  17),  tho 
experimental  and  analytical  determination  of  the  mechanics  of  rarefied 
Rases  is  being  conti.  uaily  advanced.  Consequently,  it  is  expected  that 
satellite  aerodynamics  in  the  near -Earth  regime  (say  less  than  hOG  miles 
altitude)  will  warrant  increased  attention  as  more  sophisticated  vehicles 
are  inherently  utilised. 

The  application  of  rarefied  gas  kinetics;  viz.,  free-molecule  aero¬ 
dynamics,  has  beer,  until  recently,  hopelessly  complicated  for  all  but  a 
few  academic  problems.  These  classical  applications  (e.g.,  Refs.  8,  17, 
18),  however,  have  provided  foundations  which  have  led  to  the  development 
of  approximate  free-moleculc  methods  for  estimating  aerodynamic  forces 
and  moments  on  convex  bodies.  Cne  notable,  but  restrictive,  approach  is 
the  rNewtonian -Diffuse"  method  of  Gustafson  (Ref.  19,  20)  which  implies 
that  the  random  molecular  motion  is  insignificant  when  compared  with  the 
vehicle* s  motion.  Another  ^.ore  general  free-molecule  theory  (Ref.  6) 
was  developed  which  considerably  simplified  the  force  and  moment  calcula¬ 
tions  vitnout  seriously  compromising  the  accurary  of  the  results.  This 
method,  which  is  outlined  and  used  herein,  only  requires  that  the  average 
normal  component  of  the  molecular  speed  ratio  be  greater  than  unity  over 
the  portion  of  the  body  Surface  exposed  to  the  free-strean.  For  roost 
satellite  configurations,  however,  this  requirement  is  compatible  with 
the  expected  conditions.  Moreover,  there  are  no  assumptions  as  to  the 
nature  cf  the  molecule-surface  interaction  process,  other  than  the  inter¬ 
action  may  be  described  by  the  usual  momentum  exchange  coefficients, 

<T  and  <r'  (Ref.  5);  and  in  fact,  these  reflection  coefficients,  as 
well  as  the  satellite  surface  temperature,  may  be  regarded  as  variables 
over  the  exposed  body  surface. 

As  in  the  case  of  any  aerodynamic  configuration,  the  aeroctynamic 
design  of  a  satellite  must  consider  the  external  shape  and  surface  condi¬ 
tion  of  the  vehicle  as  well  as  the  internal  mass  di  *t*ibvtion,  or  inertia 
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characteristics.  This,  of  course,  means  that  certain  design  compromises 
are  inevitable  among  the  various  aevices  sbcard  the  vehicle;  e  £.,  solar 
"paddles”  may  not  oe  compatible  with  efficient  aerodynamic  design.  Never¬ 
theless,  feasibility  studies  of  satellite  aerodynamic  stabilization  must 
oe  kept  elesenlary  anc  uncluttered.  It  is  for  this  reason  that  the 
discussion  in  this  report  is  primarily  concerned  with  simple,  aerodynamic 
configurations.  These  basic  designs,  however,  are  certainly  no  more 
idealized  than  the  classic  "dumbbell"  satellite  used  in  gravity  gradient 
studies. 

Although  this  study  is  primarily  concerned  with  the  pitch  stability 
arid  orientation  of  specific  satellite  configurations,  aej odynamic  force 
and  moment  equations  are  presented  (Ref3.  b.  6)  which  may  be  utilized  for 
calculating  aerodynamic  coefficients,  both  static  and  dynamic,  of  other 
convex  body  shapes  not  treated  herein.  In  addition,  the  basic  free- 
moiecule  force  ana  moment  equations  are  amenable  to  the  calculation  of  the 
effects  of  solar  pressure,  as  illust: ated  in  Appendix  A  of  this  report. 
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1.  REVIEW  of  FdEE-HOlSCULB  AKlOSTHtfL'CS 


1.1  Introduction 

Tho  purpose  of  this  section  is  to  review,  without  elaboration,  the 
nomenclature  and  equations  which  are  frequently  referenced  throughout 
this  volume.  Most  of  this  material  has  been  presented  extensively  in  the 
open  literature,  and  consequently,  original  sources  and/or  more  detailed 
analyses  will  be  cited  whenever  possible. 

1.2  Shear  and  Pressure  Coefficients 

For  a  surface  element  exposed  to  a  free- molecular  stream,  the 
impinging  molecules  do  not  strike  the  rebounding  molecules  near  the  sur¬ 
face  with  sufficient  frequency  to  vary  the  incident  momentum  and  energy. 
This  hypothesis  therefore  permits  separation  of  the  incident  momentum 
f-ux  ,rom  the  momentum  flux  due  to  the  reflected  or  re-emitted1  molecules. 
Using  the  standard  approximations  for  a  molecular  stream  m  Maxwellian 
equilioriun,  the  surface  tangential  (snear)  and  normal  (pressure) 
components  of  the  incident  momentum  flux  per  unit  area  may  be  written  as 
(dels.  1  through  5,  e.g. ) 


(2) 


where  S-ri  and  Sr  are  respectively  the  normal  and  tangential  components 
of  the  free-strean  molecular  speed  ratio 

The  incident  pressure  and  shear  stresses  determined  above  must  now 
bo  combined  with  the  contributions  of  the  reflected  molecules  to  give 
the  total  free-Kolecule  stresses.  The  conventional  (but  not  universal) 


1Because  of  the  empirical  natura  of  the  analysis  to  follow,  a  discussion 
of  the  conceptual  difference  between  "reflection"  and  "re-enission"  is 
not  warranted. 


h 


practice  .s  to  rebate  ti  e  reflected  shear  and  pressure  stresses  to  the 
incident  stresses  by  th#  empirical  reflection  coefficients  cr  ard  cr’ 
(Re^s.  1,  ?  and  5)«  i.ieae  coefficients  describe  the  areiunt  of  nonentum 
trar.3Tj.tfcea  to  the  codjy  surface  in  ter.~s  of  cepartures  : 'top  completely 
diffuse  reflection.  At  tus  extreme,  the  molecules  are  assumed  reflected 
with  KiAwellian  random  rotior.  corresponding  to  tne  surface  temperature. 

At  the  other  extreme,  the  incident  molecules  are  assumed  to  be  specularly 
reflected  with  the  normal  component  of  velocity  unchanged  m  magnitude  hut 
reversed  in  direction.  Then,  by  using  tr.ese  extremes,  the  reflection 
coefficients  are  defired  by  (Ref.  1,  e.g.) 


u  -  T  r 

\  -  Ti 


*  -  Pd 

where  the  subscript  r  refers  to  the  actual  stress  due  to  the  reflected 
molecules,  and  where  the  subscript  A  refers  to  the  value  of  the 
reflected  stress  cue  to  wholly  diffuse  reflection.  Then  when  «*=  c ' =  1 
the  reflection  is  diffuse,  and  when  O'  =  <r‘=  o  ,  the  reflection  is 
specular.  The  portion  of  the  shear  stress  due  to  diffuse  reflection  (Td)# 
however,  is  identically  cere  since  there  is  no  preferred  tangential  direc¬ 
tion  for  molecules  reflected  diffusely.  On  ,he  other  hand,  the  diffuse 
pressure  term  p.  has  a  non-zero  value  and  has  been  derived  as  (Refs.  1 
and  3) 


a  feints  'l'*-*5"!} 


(5) 


where  Tw  and  T*»  are  respects  ely  the  surface  and  fVee  stream  tempera¬ 
tures. 

By  combining  Sqs.  (l)  through  (5),  it  is  easily  shown  that  the  total 
shear  and  pressure  coefficients  are 
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By  th t-  method  of  Reference  h  or  6,  the  Gaussian  error  function  in 
the  foregoing  eouationsnay  be  rearranged  and  integrated  by  parte  to  give 
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where  the  last  tern  nay  now  be  regarded  as  a  snail  remainder  of  a  two* 
tern  approximation.  For  a  surface  elenent.  exposed  to  the  free-strean 
lieu,  substitution  of  £q.  (16)  into  Eqs.  (6)  and  (7)  yields 
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where,  for  greater  their,  unity, 
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For  a  surface  element  shielded  from  the  frer  -stream,  S*j  and  5T  may  be 
considered  negative  and  Eqs.  (6)  and  (7)  become 


cp,  .  tflf* 

a-a-1  e-s»  (13) 
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where  it  should  be  noted  that 

erf(~x^  =  ~erf(x^  • 

If  the  order  of  magnitude  of  Eq.  (11)  xs  utilized,  then  Eqs.  (9)  and  (l<>) 
are  approximately  (Ref.  6), 

Cr  =  zrJr/SjA*-  (ll|) 
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Likewise  for  the  shielded  cr  shadowed  element. 
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(16) 
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Eos.  (Hi)  through  (17)  are  based  cn  the  for ro}  restriction  that  |S#f(  —  1 
at  every  point  on  the  body.  For  most  applications,  however,  it  is  suffi¬ 
ciently  accurate  to  assume  that  the  foregoing  equations  are  valid  for 
the  gpan  on  the  body  greeter  than  unity. 

1.3  Force  and  Mcoent  Equations 

Force  and  pitching  moment  integrals  may  now  be  formulated  using  the 
stress  coefficients  obtained  above.  Referring  to  Fig.  1,  it  is  not 
difficult  to  show  that  the  drag,  lift  and  pitching  moment  coefficients 
for  a  surface  element  are  (Ref.  ii) 
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where  Si  and  o  are  the  reference  area  and  length  respectively.  For 
an  arbitrary  body  of  revolution,  it  ia  elenentary  to  show  that 
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where  f*  -  **7 Jx  and  the  nornal  vector  N  is  always  directed  anwtrd 
(convex  body).  The  substitution  of  ')qs.  (19)  into  Eqs.  (lb)  yaelcs  the 
fdU'ving  general  expressions  for  tlu  force  and  Kouent  coefficients  for 
a  bocy  of  revolution: 
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These  expressions  ceco.-ne,  upon  substitution  of  the  approximate  stress 
coefficients. 
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Cn. 


+  2.<T[|+  r‘x 


«sc<  +  (x-x£m.\sin<<]| 


(21) 

(cent.) 


vr.ere  the  integration  now  annlies  only  on  the  portion  of  the  body  exposed 
to  the  incident  molecular  stream.  Referring  to  Fig.  1  ar-a  2qs.  (19),  it 
may  be  seer,  that  for  the  steady  state  case, 
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In  a  similar  manner,  the  incremental  velocity  distribution  can  be  found 
for  a  body  patching  about  its  center  of  mass:  i.e.. 


rr'+  A  -  X „ 


(23) 


where  %  is  the  pitch  rate.  The  effective  normal  component  of  the 
molecular  speed  ratio  is  then 


where  j 2  V*>  is  the  c  cnventional  non-dimensional  pitch  rate . 

The  approximate  static  force  and  moment  coefficients  for  an  arbitrary 
body  of  revolution  are  found  by  substituting  Eq.  (22)  into  Eqs.  (21).  The 
result  is 
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where  X  is  the  body  length  and  where 

=  (~r '/iom*)  ‘fw  r'  1  |-tav.4\ 

-  IT  -for  r‘  S  . 


*For  integratio 1  purposes,  X  is  nov  being  defined  as  positive  aft  of  the 
nose. 


An  aerodynamic  parameter  frequently  utilized  is  the  static  pitchirg 
m onent  coefficient  slope  at  zero  angle  of  attacK.  For  bodies  of  revolu¬ 
tion  with  r'  >  o  >  the  derivative  of  the  third  of  Eqs.  (?5)  yieius 

“  5nc)  ]  +  r.*  2S.^L  f  H 

V°  0  (26) 

+iT  fr+2r'(x->w)j|  rdy 

where  O',  0*  *  and  Tw  are  taken  to  be  independent  of  the  azimuth  angle 
$  .  The  dating  in  pitex.  derivative  Cm-,  may  be  found  by  differenti¬ 
ating  the  third  of  £qs.  (21)  with  respect  to  ac/ZVao  while  using 
Eq.  (2i*);  i.e., 
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The  foregoing  integrals,  Eqs.  (16),  (19/,  (25/  through  (28),  wil.  be 
useful  in  the  <irag  and  moment  analytes  to  follow.  When  possible,  the 
oasic  Eqs.  <18)  and  (20)  will  be  solved  using  the  exact  shear  and 
pressure  coefficients  Eqs.  (6)  and  (7).  For  most  applications,  however, 
the  approximate  Eqs.  (25)  through  (27)  will  satisfy  the  accuracy 
requirements. 
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DRAG  ANALYSIS 


2.1  Introduction 

In  the  expressions  for  satellite  lifetime  and  atmosphere  influence 
on  orbital  elements,  the  well-known  ballistic  coefficient.  C&Sir/Zm- 
inherently  appears  .  The  problem  of  evaluating  the  coefficient  may  be 
phrased  two  ways:  if  the  drag  coefficient,  Co  ,  is  based  on  a  fixed 
reference  area,  S*  ,  then  Co  is  a  strong  function  of  the  attitude  of 
the  satellite;  alternately,  if  an  average  Q>  is  assumed  for  all 
satellite  attitudes,  then  the  proper  reference  area  becomes  variable. 

It  is  this  latter  approach  which  is  more  coraaonly  taken.  That  is,  let 
Sr,  equal  the  projected  frontal  area  of  the  satellite  (normal  to  free- 
strearn),  and  then  assume  a  constant  Cd  which  is  based  on  this  frontal 
area.  The  discussion  in  this  section  is  concerned  with  the  verification 
or  nullification  of  this  assumption.  Three  practical  configurations  - 
cylinder,  cone  and  sphere  -  are  considered  along  with  the  "academic" 
flat  plate. 

All  free-molecule  drag  coefficients  are  calculated  using  the  exact 
methods  when  possible.  The  approximate  method  of  Ref.  6,  which  is 
outlined  in  Section  1,  is  used  for  comparisons  or  when  the  exact  expres¬ 
sions  are  not  integrable.  The  reflection  coefficients,  <T  and  (T1  ,  as 
well  as  the  surface  temperature  Tw  are  assumed  constant  over  the  body 
surface.  If  the  free  stream  temperature  T-o  is  taker,  from  the  ARX  1959 
Model  Atmosphere  (Ref.  7),  and  if  the  relative  velocity  is  assumed  to  be 
approximately  the  inertial  velocity,  then  L.e  variation  of  the  speed 
ratio  with  altitude  for  nearly  circular  orbits  may  be  calculated. 
This  variation  of  ,  along  with  a  curve  of  (l/Sm •)  t7W/r«»  for 
T*  =  550°R  ,  are  shown  in  Fig.  2  and  are  used  in  the  drag  analyses  to 
follow. 

2.2  Cylinder 

The  cylinder  drag  coefficient  must  be  divided  into  two  parts:  that 
due  to  the  blunt  ends  and  that  due  to  the  cylindrical  surface.  The 
corresponding  frontal  areas  are 


Sw  =  UR’-MS* 


-  TIZ*' jjtos “I  + 


See  Ippendix  0  of  Volume  I  of  this  series  of  reports. 
^Gustafson  (Ret.  ?0)  also  asrunes  "  5$0°R. 


yjicre  all  starred  ayaibols  will  now  refer  to  the  frontal  area  syst-en.  The 
free-molecule  drag  coefficient  of  the  coabined  f)at  nose  and  base  is  the 
sane  as  a  flat  plate  at  angle  of  attack.  This  result  nay  be  found  as 
(Ref.  2  or  h)  _ 
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where  the  third  of  Eqs.  (29)  is  used  for  the  reference  area.  Tbs  free 
molecule  drag  coefficient  of  the  cylindrical  portion  may  be  formulated  in 
a  similar  manner  except  that  the  integrations  are  more  tedjous.  The  exact 
expression  has  been  obtained  by  Talbot  (Ref.  8).  The  result  is. 
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where  |  =:  5t*Sw  cv  -  and  where 
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If  diffuse  molecular  reflection  is  assumed,  then  ^=0"’=  |.o 
and  Eq s.  (30)  and  (31)  may  be  combined  and  simplified  to  yield  the  total 
drag  coefficient. 
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How  for  near-Earth  mbits,  Swo  is  always  in  excess  of  5,0;  therefore, 
e-5.3--^c  and  e.rf(S^_^|.0  • 
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which  is  sillily  greater-  than  2  for  snail  2.(0  ,  For  large  JtlQ  , 
novever ,  the  cylinder  "skin  friction"  tern  4(lf0)/fir  S<x>  nay  be 
doninant.  It  should  be  noted  th3o  the  approximate  theory  of  Ref.  6  is 
not  applicable  for  this  ©t-o  case  since  *>*>  on  the  cylindrical  portion 
is  zero.  The  first  three  terns  of  Eo.  (35),  however,  are  due  to  the 
blunt  nose  and  would  be  obtained  frcra  the  approximate  theory. 

At  another  extreme,  ci  -Ti/t.  ,  the  reference  area  becomes  J^D  and 
Eq.  (33)  reduces  to 
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This  equation  nay  also  be  simplified  for  large  since  it  can  be 

sncwn  that  ¥($+)  -»  2 ;  i.e.. 
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which  is  seen  zo  be  near  2  for  all  reasonable  AlO  and  •  Conse¬ 
quently,  for  large  ilD  *  G?  at  angles  of  attack  of  zero  and  --ay 
be  considerably  different. 


2q.  (33)  may  be  solved  for  other  values  cf  angle  of  attack  and  fine¬ 
ness  . atio.  Several  combinations,  along  with  the  approximate  Eqs.  (35) 
and  (37),  are  «:icwn  in  F5g.  3  for  C *  versus  .  It  is  interesting 
to  note  from  fig.  3  that  the  minimum  r ?  of  a  cylinder  occurs  for  some 
intermediate  angle  of  attack. 
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2-3  Cone 


Unfortunately,  the  exact  free-noleculc  drag  coefficient  has;  not 
been  explicitly  determined  for  the  case  of  a  right  circular  cone  at  angle 
of  attack.  Nevertheless,  numerical  solutions  have  been  obtained  for 
inclined  cones  and  are  presented  by  stalder  and  Zurick  (Ref.  9)  for 
completely  diffuse  and  completely  specular  reflection.  For  zero  angle  of 
attack,  however,  the  first  of  Eqs.  (20)  is  easily  integrated  to  give  the 
exact  drag  coeflicient  as 
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where  S  is  the  cone  sem vertex  angle  and  where  tne  base  contribution  is 
neglected.  For  diffuse  molecular  reflection,  T--  <T'  c  |  o  and  Eq.  (38) 
reduces  to 


As  becomes  larger  than  unity,  the  exponential  and  error 
function  approach  the  limits  zero  and  one  respectively,  as  previously 
noted,  and 


ft  S  iL|]t  +  2  +  -A-  UO) 

* -o  i  ico 


which  is  the  same  as  the  solution  of  the  first  of  Eqs.  (25),  for  O  . 
For  very  thin  cones,  Eq.  (39)  predicts  free -molecule  drag  coefficients 
considerably  higher  than  calculated  by  Eq.  (liO),  as  may  be  seen  in  Fig.  1*. 
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for  cones  at  angle  or  attack,  the  approximate  theory  outlined  in 
Section  1-3  (or  Ref.  6)  may  be  cc: veniently  utilized  if  the  average  SfO 
on  the  cone  is  greater  than  unity .  This  restriction  may  be  expressed  as 


s„>  _ 1___ 

n(<)  C-oSoO  S| Y\  & 
where  £  and 

K(*\  a  -ki  ^''(-<)  +  fTT^  L^-i0] 

n  dii) 

=  i  o 

pit)  =  [<ilo] 

=  I O  f  £  k  I  o  j  •  (ii2j 


Since  the  frontal  area  of  the  inclined  cone  must  also  be  known,  it  is  not 
difficult  to  show  that 


vnere  R.  13  the  radius  of  the  base.  The  first  of  Eqs.  (25)  may  be  inte¬ 
grated  and  combined  witn  do,s.  (U3)  to  give 


Q>  =  -  +  jljjLA—  inaito 

w(*f  )  s»  h«  v>«) 

-fizl'  +  zr 
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triier*' 


Hi) *  +  j( ' I  < 1 !  !<s  1  ‘l 

«i-*5V,  [<>->«]  “ 
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=  i  + 
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Eq.  (Itl»)  is  also  presented  in  Refs,  lx  and  6  for  the  case  oi  =  TrR?-. 
For  diffuse  reflection,  Eq.  (!x!x)  sucplifies  to, 


c;s  i;1%  «<“*#>* 2  *<h 


(1*7) 


For  cl  -  TT/X  , 


$!*S  C<JS<< 


■»«) 


and  Eq.  (Ij7)  becomes. 


cSL,  s  +  n-L 

4S»  It., 


(u8) 


Eqs.  (39),  (1x0),  (la 7)  and  (U8)  are  plotted  in  Fig.  lx  for  several  angle 
of  attack  and  cone  angle  combinations. 


*i’he  last  terra  of  the  cone  Co  equation  in  Ref.  6  should  be: 

(lz^i+2rT)cc5«  k(<) 
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2.2i  Sphere 

The  sphere  is  an  extremely  simple  case  to  analyze  since  the  frontal 
area  remains  constant  regardless  of  angle  of  attack.  The  exact  free- 
molecuie  drag  coefficient  for  a  sphere  is  vritten  as  (Ref.  1,  2  or  U) 


g  -  (z -r'+s)<X$  + 

3 


( h9 ) 


wVore 


= 


z^  +  ! 

2(5  »  *■  ^ 


4-C  +4C  -  I 

4S  * 


eK(S„')  . 


(50) 


Eq.  (h?),  for  diffuse  reflection  is  obviously 


Co*  -  ZQfa  4-  Ijf  0~  . 


(51) 


"he  approximate  drag  coefficient  of  a  complete  sphere  may  be  obtained  by 
integrating  the  first  of  Eqs.  (25)  over  the  forward  hemisphere;  i.e.. 


(2 


S  2+  Jl-  + 
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Z*TT  fjr. 

5S«  I  TX 


(52) 


where  the  shadow  area  contributes  very  little  to  the  total  drag.  The 
difference  between  these  last  two  equations  may  be  seen  in  Fig.  f>.  Her . 
again,  the  drag  coefficient  approaches  the  "Newtonian-diffuse"  value  of  2. 

2.5  Flat  Plate 


The  infinitely  thin  flat  plate  is  difficult  to  analyze  from  the 
frontal  area  point- cf -view  since  at  zero  angle  of  attack  the  area  is  zero 
while  the  drag  is  finite;  i.e., 


Co*  °° 


~W — "T 


1" — ▼  a  ”  i  3  j*  u  .  “5  r~5i  ^  Q  o  “^r 
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It  is  interesting,  however,  to  calculate  for  other  angles  of  attack 

for  coup  arisen  with  the  three-dimensional  configurntiom,  Ict  arbitrary 
values  of  tne  reflection  coefficients,  the  exact  Cg  may  be  expressed  as 


C?  •=  (T  fusing 
S- 


f£  +  )l£ _ +  z(*-f 

'  1  T»  fir  Sow^wci  'Va.ft  5 

4^(lft  *~<r'\  4-  z (2-r-f ;)410-<<Z 

Srf*  )  J 


v,hich  becomes  for  CT  ‘T 1  1 . 0 


C*  -  fir  vto«(  'Z.e-  a-  (7  t-_L\ 

5*  iT«  ••  S“' 


and  for  (««  larger  than  unity. 


£*  ~  fv  tloot  I1!*  ^  2. (55) 

5«  I T*»  5^ 

A  flat  plate  normal  to  the  free-strear.  («t  -  <3 o')  has  an  exact  diffuse 
drag  coefficient  of 

til  +  <*» 


which  approaches,  for  large  speed  ratios. 


£?  a  ffe  +  It  -h  ■ 

•('-■J/2.  s„  1T«  ’•* 


Eqs.  (Sir)  through  (57)  arc  plotted  in  Fig.  6  for  several  angles  of  attack. 
The  deviations  at  snail  angle  of  attack,  as  for  the  case  of  the  cylinder, 
are  clear rv  seen. 


2.6  Surgery  or  Configurations 

A  comparison  of  diffuse-  frontal  drag  coefficients  for  one  cylinder, 
cone,  sphere  and  flat  plate  is  shown  in  Fig,  7.  If  onlj  the  three- 
dimensional  shapes  are  considered,  it  nay  oe  seen  that  most  of  the  drag 
coefficients  lie  in  a  t$%  band  for  all  practical  molecular  speed  ratios. 

A  noteworthy  exception  is  the  cylinder  rear  2ero  angle  of  attack.  Then 
it  nay  be  concluded  that  the  practice  of  assuming  a  constant  drag 
coefficient  per  unit  frontal  area  for  a  satellite  at  all  attitudes  is 
not,  in  general,  correct,  although  for  seme  shapes  it  may  be  an  acceptable 
approximation.  The  problem  of  estimating  the  ballistic  coefficients, 
therefore,  appears  to  deserve  additional  study. 
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?.  STATIC  STABILITY  ANALYSIS 


3.1  General  Considerations 


4s  noted  in  Volume  I  of  this  series  of  reports,  the  aerodynamic 
contribution  on  the  equations  of  satellite  pitching  motion  may  be  conven¬ 
tionally  represented  by  the  leading  terms  of  a  Fourier  series.  For  small 
angular  deviations  of  the  body  from  equilibrium,  it  may  therefore  be 
shown  that  the  non-dimensional  aerodynamic  torque  ior  pitching  moment 
coefficient)  may  be  written  approximately  as 

U  +Cn(-^\  (58) 


where,  in  this  analysis, 


v° 


j  fit  \ 

t^JJ 


d~Q 


(59) 


The  discussions  in  this  section  are  concerned  only  with  the  static 
pitching  moments  (i.e.,  with  O  )  and  the  initial  pitching  moment 
slope  Cm«t  •  The  purpose  of  the  first  of  these  investigations,  namely 
Crw  versus  a  wide  range  of  angle  of  attack,  is  to  examine  the  validity 
of  the  linear  Crw  approximation  for  three  simple  configurations:  t\  » 
cylinder,  sphere  and  cone.  The  analysis  is  subsequently  extended  to 
calculations  of  Cm<<  for  more  complex  configurations:  power  law  odies, 
truncated  cones,  cone-cylinders  and  spherp-capped  cylinders.  Supplemental 
equations  and  graphs  are  presented  for  the  purpo  of  examining  miscel¬ 
laneous  effects  of  surface  temperature  and  reflection  coefficient  varia¬ 
tions.  The  final  portion  of  this  section  is  concerned  with  the  optimiza¬ 
tion  of  satellite  geometry  for  maximum  aerodynamic  stability.  The  possible 
au mentation  effecus  of  gravitational  torques,  however,  are  not  ignored. 


30 


3.2  Aerocynamc  Pitching  Me  rent 
3-2.1  Introduction 


With  proper  satellite  design,  numerical  results  of  the  motion 
analysis  (Volume  I)  have  sr.ovn  that  the  small  angle  of  attacK  approxima¬ 
tion  is  generally  justified.  In  order  t«  avoid  any  unforeseen  complica¬ 
tions  ,  however,  it  is  advisable  to  examine  the  pitching  characteristics 
beyond  the  small  angle  range.  Acccraiiigly,  it  is  the  purpose  of  this 
section  to  present  curves  of  pitching  moment  coeflicient  versus  angle  of 
attack  for  the  oasic  oody  shapes  *v  cylinder,  sphere  and  cone).  Initial 
pitching  nceicnt  slopes  for  tnese  shapes  are  also  presented  along  with  the 
slopes  for  the  more  sophisticated  shapes  (power  la«  bodies,  truncated 
cones,  hem: sphere -cylinders  and  cone-cylinders). 

Ail  coefficients  are  referenced  to  the  base  diameter  O  and  base 
area  it  dx/4  • 

3.2.2  Pitching  Moment  of  Baric  Body  Shapes 
3.2.2. 1  Cylinder 


The  pitching  moment  coefficient  of  a  right  circular  cylinder  may 
be  calculated  using  the  exact  free  molecule  flow  e  ,  tiens.  For  a 
cylinder  without  end  effects,  the  last  of  Eqs.  (20;  integrates  to  give 
(Ref.  1) 
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where  f  and  ?*(?)  is  defined  by  the  first  of  Sqs.  (32). 

The  nose  contribution  may  be  taken  into  account  by  considering  a  flat- 
olate  normal  to  the  stress.  The  exact  Moment  contribution  of  tnis  surface 
is  due  only  to  the  tangential  stresses  and  is  given  by 

C"We=  fjff-  T *"*['  +  •  <61J 


It  is  noted  that  the  nose  contribution  is  positive  and  hence  destabilizing. 
The  base  end  effects,  which  are  snail,  are  not  included  in  this  analysis 
since  the  base  will  be  assumed  v>id  fer  design  purposes  (see  Section  3.ii). 
CoRoming  Eqs.  (60)  and  (6l),  vhe  total  static  pitching  nonent  coefficient 
is 


c, 


O-'fiT 


+  sm2c(  ^er-f(s<Dtos^  -  lj"^ 


(62) 


For  diffuse  molecular  reflection,  <T- <T  ‘  -  1  O  >  and  £q.  (62)  becomes 


Eq.  (63)  as  shewn  an  Fags.  13,  12  and  13  foi  cylinder  fineness  ratios  of 
i,  ?  ard  5  respectively  and  for  Tv  -T-o  .  The  temperatures  are  equated 
to  simplify  the  presentations.  The  effect  of  Tv  £T«®  >  however,  is 
discussed  an  Section  3. 3. 2.1. 

It  may  be  noticed  that  the  last  term  in  Eq.  (63)  is  quite  small 
when  is  in  excess  of  about  2.0.  Then  for  small  angles  of 

attack  this  term  may  be  neglected.  The  pitching  moment,  slope  at  zero 
angle  of  attack  may  be  foritCG  by  differ  entiatuig  Eq.  (63)  with  respect  to 
oL  and  passing  to  the  limit  <*-->0  .  Accordingly,  one  finds. 


where  the  last  term  in  Eq.  (63)  is  ignored  and  where 


From  Eq.  (61*)  the  cylinder  is  seen  to  be  aerodynamically  stable  for  all 
centers  of  mass  forward  of  the  geometric  center.  The  results  of  Eq.  (61* ) 
will  be  utilised  in  the  cylinder  design  analysis  of  Section  3.1*. 2.1. 

3. 2. 2. 2  Sphere 


Figure  9»  Sphere  Nomenclature 
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The  pitching  moment  coefficient  of  a  sprere  is  easily  calculable 
since  the  result  in  t  aerodynamic  force  (drag)  acts  through  the  sphere 
center.  That  is, 

Cwv.  “  “  ~  1|Q>  5iv\o<  .  (66) 

The  exact  drag  coefficient  for  arbitrary  molecular  reflection  is  given 
by  Sq.  (U9)  as 

Co  =  (z-  r'+r)^)  +  zr'K  (67) 

2>  I 

where  ($<»')  ‘‘s  defined  by  Eq.  (50).  Then  for  diffuse  leflection, 

Eq.  (66)  becomes  upor  ihe  subrtitution  of  Eqs.  (50)  and  (67) 


which  varies  linearly  with  center  of  mass  location.  The  normalized 
version  of  Eq.  (66)  is  shown  in  Fig.  lh  for  Tw  -To*  -  Temperature 
effects  are  considered  in  Section  3.3. 2.2. 

The  sphere  pitching  moment  slope  at  zero  angle  of  attack  is 


c 


v 


\Zj£+l_ 

UfiT« 


-C 

t  4 


(69) 


which,  like  the  cylinder,  is  2ero  for  coincident  centers  of  volume  and  •* 
mss. 
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Figure  10.  Gone  Nomenclature 


Unlike  the  cylinder  and  the  sphere,  there  appears  to  be  no 
straight-fcn- arc  I'.ethod  for  obtaining  the  exact  free  molecule  pitching 
moment  of  the  cone  at  angle  of  attack.  Using  the  last  of  Eqs.  (25), 
her ever,  approximate  pressure  ana  shear  coefficients  may  be  integrated 
with  lit+le  difficulty.  The  x-esulting  static  pitching  moment  coefficient 
may  be  shorn  to  be  (Ref.  h) 


r  f* 
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+  —  Jj£l  “if ,  ,  -  1 
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(70) 
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i.he-c  tie  parameter  <  .s  defined  as  i*ni/ and 
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The  foregoing  functions  are  the  result  of  integrating  tne  pressure  and 
shear  over  the  surface  area  exposed  to  the  incident  molecular  stream. 
The  contribution  of  the  "shadov"  is  neglected  in  the  approximate  theory, 
I.  u  therefore  clear  that  vhen  the  angle  of  attack  is  less  than  the 
cone  angle,  the  entire  conical  surface  "sees"  the  free  stream  and  the 
functions  assume  integer  values.  Consequently,  for  diffuse  reflection, 
■A •  I/O)  is  expressible  in  two  parts: 
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+sm«  ios«  jft<)  - 

2qs.  (72)  are  shown  plotted  in  Figs.  1$,  16,  and  17  for  cone  angles  of 
l6c»  30°  and  ltp°  respectively  and  for  Tv,  -  .  It  nay  be  noticed  in 

these  figures  that  there  is  a  slight  discontinuity  in  the  pitching  moment 
slope  at  the  point*  where  <*  equals  £  .  These  discontinuities  are  due 
to  the  sudden  appearance  oi  a  shadow  as  the  angle  of  attach  increases 
beyond  the  eone  angle.  If  exact  fre^-raoiecule  theory  were  available  for 
these  conditions,  hovev  ,r,  the  shadow  ax.^a  would  be  properly  taken  into 
account,  and  the  discontinuity,  if  present  at  all,  would  be  »noothed. 

The  static  pitching  norisnt  derivative  at  zero  angle  of  attack 
may  be  found  from  Eq.  (?6)  or  by  differentiating  the  first  of  Eqs.  (72) 
with  respset  to  o(.  and  taking  the  limit  as  o i  approaches  zero.  The 
result  is. 
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which  is  applicable  for  So-  sm  £  greater  than  unity. 
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3.2.3  Pitching  Moment  Slopes  of  Cxiplex  Body  Shapes 
3.2. 3.2  Power  Law  Body 


Figure  18.  Power  Law  Body  Nomenclature 


As  the  result  of  the  cylinder,  sphere  and  cone  stability 
comparison  of  the  previous  section,  it  seems  worthwhile  to  continue  the 
analysis  for  somewhat  more  practical  shapes.  The  power  law  body  is  such 
a  shape,  and  hence,  is  described  below. 

A  power  bod}  of  r*  vclvtc'W'  is  described  by  the  well-known 

formula, 


where  £  and  H  are  the  maximum  body  radius  and  1  »th  **espectively,  and 
where  O  <  v\  ^-1  .  Now  the  slops  of  the  pitching  moment  coefficient. 

Cm*  .  at  zero  angle  of  attacK  for  a  body  of  revolution  with  r’  ^  O 
may  be  determined  from  Eq.  (26)  using  the  approximate  free-mo3ecule  theory. 
The  result  of  substituting  Eq.  (7h)  and  its  derivatives  in  Eq.  (26)  is 
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Except  for  a  few  special  cases,  the  first  two  integrals  in  Eq.  (75)  have 
not  been  evaluated  (at  presen:)  for  a  general  value  of  n  .  For  diffuse 
molecular  reflection,  howeve’*,  (T-(r‘  =  1.0  and  the  coefficient  of  the 
first  integral  is  identically  zero.  Also,  sai-ple  calculations  have 
shown,  using  Fig.  2,  the*.  tne  temperature-dependent  term  is  less  than 
of  the  dominant,  last  terra.  Therefore,  if  diffuse  reflection  as  assumed, 
and  if  the  second  term  is  neglected,  Eq.  (75)  becomes 


This  integral  rasy  be  evaluated  exactly,  giving  the  result, 

which  is  independent  of  R  .  Consequently,  for  diffuse  reflection,  the 
only  effect  of  the  exponent  on  pitching  moment  slope  is  due  to  the  small 
contribution  of  the  reflected  pressure  coefficient,  as  treated  in  Sec¬ 
tion  1.  Also,  it  may  be  seen  that  Eq.  (76)  agrees  with  the  leading  tern 
of  Eq.  (7;)  for  the  cone  =  • 


2.2. 3.2  Truncated  Cone 

The  stability  analysis  of  the  more  complex  geometrical  shapes 
continues  herein  for  the  case  of  a  cone  frustum.  The  calculations  are 
considerably  more  complicated  for  this  configuration  sir.ce  the  degree  of 
tn  cation  becomes  an  additional  variable.  The  geometry  and  nomenclature 
used  for  the  subsequent  calculations  are  shown  below  in  Figure  19. 
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Figure  19.  Truncated  Cone  Nomenclature 


The  free -molecule  pitching  moment  slope  of  this  shape  may  be 
found  by  combining  the  separate  contributions  of  the  conical  surface  and 
the  flat  nose.  Aitnough  the  exact  free -molecule  flow  theory  may  be  us  d 
for  the  nose  portion,  it  is  sufficient  (and  convenient)  to  calculate 
Cm«  at  2ero  angle  of  attack  using  the  restrict'  “*s  tha*  5®  >6 
and  >  1  .  Then  for  diffuse  molecular  reflection,  Eq.  (2o) 

shows  that  Cm*  of  the  conical  portion  ( dr/Jx  /  O ^  is 


)i7 


and  correspondingly,  the  nose  contribution  is 


where  £  is  the  ratio  of  the  nose  diameter  to  the  base  diameter.  The 
total  static  pitching  ronent  slope  is  the  sure  of  Eqs.  (77)  and  (76). 
That  is. 


It  is  not  difficult  to  show  for  reasonable  configurations  with 

)  2.  that  t-he  temperature-dependent  tars,  in  £q.  (79)  is 
considerably  smaller  than  the  last  term.  Although  the  difference  will  be 
shown  later,  it  is  now  convenient  to  neglect  the  temperature  term  and 
write  Eq.  ( f$)  as 


(So) 


which  is  the  sane  equation  used  in  the  previous  section  for  power  law 
bodies.  As  Eq.  (80)  is  presented  above,  Cm*  is  independent  of  the 
diameter  ratio  0  . 
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3.2.3.'*  Henisphcre-Cy  lander 

The  geometry  of  the  herdsph ere -cylinder  is  shown  in  Fig.  20. 


Figure  20.  Hemisphere-Cylinder  Nomenclature 


The  pitching  moment  slop?  of  a  herd  sphere -cylinder  may  be  found 
by  adding  the  spuei^cal  cap  contribution  to  the  cylinder  contribution, 
both  of  which  have  been  discussed  previously.  As  approaches  zero  the 
value  of  Cm*.  for  the  hemisphere  is  essentially  the  same  is  the  value 
for  the  complete  sphere  since  the  contribution  of  the  aft  portion  of  the 
sphere  is  negligible.  Then  from  Eq.  (69) 


where  Q(S*)  is  defined  by  Eq.  ($0)  and  where  diffuse  reflection  4s 
assumed.  Likewise,  Cm*  for  the  cylinder,  without  end  effects,  auy  be 
found  by  differentiating  Eq.  (60)  with  respect  to  oL  and  passing  t 
the  limit.  The  result  is,  for  <T - T'  =  1.0  > 


wnere  «.«’  x*mc.  are  measured  from  the  origin  of  the  cylinder. 
Frcra  Fig.  ix.  is  seen  that 
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X  =  5C  +  r>/z 
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Chess  relations  naj’  be  substituted  in  3q.  (82)  to  provide  a  cc“r*.on 
reference  basis;  i.e., 

r  =  _/X-i.'5ilU  JL  \(i'U  1  +i.  -2iM  +  l?  («) 

uv,L  (0  -2 jjf?U,  k<Jz  jLd  i  d  J  y 

where  £/t>  is  greater  than  ^2  The  combination  of  £qs.  (81)  and  (83) 

gives  the  total  pitching  moment  slope  a  hemsphci  e-cylinder  with 
diffuse  reflection: 
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3*2. 3.U  Cone-Cylinder 


H-H 


The  scatic  pitching  moment  slope  of  the  cone-cylinder  configura¬ 
tion  is  easily  obtained  since  both  portions  have  beer,  previously  calculated. 
For  diffuse  reflection,  tne  approxi-oate  free- molecule  ilov  theory  gives 
for  the  conical  portion. 


r  ~  /-(fir  ImI  (*!  -  3'%s‘ 

Ok,  s  - — S. - t  \-2Mis: 


(86) 


which  is  the  same  as  Sq.  (?3)  except  for  a  slight  rearrangement.  The 
C»no(  due  to  the  cylindrical  portion  may  be  written  using  Eq.  (82)  and 
rig.  21,  noting  that 


i  =  1.  +  (tW  S 
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where  and  are  measured  aft  fro^  the  cylinder  shoulder.  Thus, 

Eq.  (&2)  beccnes 


where  ■iis  is  greater  than  Vz  .  The  total  ccr.e-cylir.der  is  the 

sun  of  Sqs.  (86)  and  (f>7): 
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3.'-  Miscellaneous  Effects  on  Stability 

3.3.1  Introduction 

As  stated  several  tines  throughout  this  report,  the  complexity  of 
free-molecule  force  and  moment  calculations  is  reduced  if  the  surface 
temperature ,  Tw  ,  and  the  tangential  and  normal  momentum  coefficients, 
CT  and  XT*  ,  are  assumed  constant  for  a  given  configuration.  Although 
those  assumptions  are  reasonable  for  most  engineering  applications,  it  as 
desirable  to  investigate  the  magnitude  of  stability  changes  caused  by 
surface  temperature  and  reflection  coefficient  variations.  In  reality, 
the  reflection  coefficients  may  vary  with  surface  temperature  as  well  as 
being  dependent  on  surface  material  and  irregularities.  For  the  subse¬ 
quent  calculations,  however,  it  will  be  assumed  that  (T  and  are 
independent  of  Tw  .  Moreover,  icr  tne  temperature  investigations,  we 
cocnonly  accepted  diffuse  values  of  tfvsfl*'  =  j.o  are  assigned,  while  for 
the  reflection  coefficient  investigations,  Tw  is  as  cun- *1  to  be  550°R 
and  Fig.  2  is  assumed  applicable. 

3.3.2  Effect  of  Uniform  Temperature  Level  on  Stability 

In  Section  3.2.2  equations  for  the  static  pitching  moment 
coefficients  were  presented  for  the  basic  body  shapes:  the  cylinder, 
sphere,  «u.d  cone.  The  accompanying  figures,  however,  were  drawn  for  the 
special  case  of  Tv.  -  To,  (for  simplicity  only).  It  is  therefore  the 
purple  wf  this  section  to  present  equations  and  graphs  for  the  incre¬ 
mental  pitching  moment  coefficient  due  to  Tv,  ^  T»  .  Since  the  basic 
equations  have  been  discussed  previously,  additional  elaboration  is  held 
to  a  minimum. 

3. 3.2.1  cylinder 


?."Ou  Eq.  (63)  t  it  is  easily  shown  that 


Eq.  (5K»  is  plotted  in  Fig.  22  versus  angle  of  attack  for  several  tempera¬ 
ture  ratios.  It  may  be  noticed  that  for  large  cylinder  fineness  ratios, 
the  temperature  ^ffe^t  may  be  rests  significant.  Also,  it  should  be  pointed 
on*  that  the  lat  nose  contribvtion  to  (shear)  is  independent  of 
temperature. 
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3. 3.2.2  Sphere 

The  effect  of  T*  £  T<©  cn  the  sphere  &*.  is  shown  to  be, 
fron  Sq.  (68) 


which  is  plotted  in  Fig.  ?3.  Unlike  the  cylinder,  temperature  only 
slightly  affects  the  sphere  pitching  reorient. 

3. 3. 2. 3  Cone 

As  indicated  previously,  the  cone  at  angle  of  attack  must  be 
treated  in  two  parts;  namely,  l«ao  tiu^x  'wits  cone  angle  £  ,  and  o ( 
greater  than  S  •  Then,  fro*  Eq.  (72), 


where  is  defined  by  the  secend  of  Eqs.  (71).  Eqs.  (92)  are 

shewn  ir.  Fig.  2h  for  cone  angles  of  15°,  30°  and  h5°.  Here  again,  as  for 
other  bodies  of  revolution  with  non-ze^o  slopes,  the  teaperature  effect  on 
stability  is  snail  as  long  as  the  fineness  ratio  is  not  large.  This 
snail  increment  nay  be  illustrated  by  plotting  the  cone  total  pitening 
mcment  slope,  Cm*  ,  versus  altitude  for  various  levels  of  surface 
temperatures.  Such  a  plot  is  presented  in  Fig.  25  using  Eq.  (73)  with  i 
cone  angle  S  =  cof (E  (an  optinum  value  as  shown  in  Section  3.H.2.3, , 

A  circular  orbit  is  assuoed  so  that  "X*  ,  and  hence  ,  may  oe 
written  as  functions  of  altitude  using  the  1959  ARI*,  toospnere  (Ref*  7). 
It  is  clearly  seen  fron  Fig.  25  that  surface  teaperature  has  a  snail 
{39?)  effect  on  the  stability  of  a  typical  cone. 


Hoi 


3.3-3  Effect  of  Asyrreetric  "'ea  per  a  tire  Varlaticr.  on  Stability 

In  order  tc  investigate  the  effects  of  an  ssy  metric  variation  of 
surface  temperature,  the  basic  pitching  reorient  integral,  Sq.  (20),  must 
be  utilized.  It  resy  be  remembered  that  for  a  completely  di  ffuse  gas- 
surface  interaction,  the  only  effect  of  temperature  is  on  the  pressure  or 
normal  momentum  component  due  to  re-eaitted  molecules  (Refs.  1,  3). 

The  ’’reflected”  pressure  coefficient  raay  be  written  frcro  Eq.  (f>) 
(for  Maxwellian  re-emission). 


where  Tv  is  regarded,  for  simplicity,  as  a  fictitious  temperature 
corresponding  to  a  surface  frem  which  the  molecules  reflect  with  complete 
thermal  accommodation. 

Since  the  shear  stress  and  the  “incident”  pressure  are  independent 
of  surface  temperature,  the  last  of  Eqs.  (20)  fer  a  body  of  revolution  at 
zero  angle  of  attach  may  be  written  as 


where  Sir  and  Z  are  the  reference  area  and  length,  and  where  the  other 
symbols  are  defined  in  Fig.  26  below. 


Figu  *e  26.  Body  cf  Revolution  at  Zero  Angle  of  Attack 


Now  consider  a  non-spinning  satellite  orbiting  at  zero  angle  of' 
attack  vath  the  son  directly  overhead.  For  a  body  of  revolution,  the 
upper  half  will  be  in  direct  sunlight  while  the  bottom  naif  is  in  a 
shadow.  If  the  satellite  is  assumed  completely  insulated  (no  internal 
heat  transfer)  and  if  free-moleuule  convpctive  h**at  transfer  is  neglected 
in  comparison  to  that  due  to  radiation,  then  the  roper  ana  lower  portions 
of  the  satellite  will  achieve  different  temperature  levels.  The  corres¬ 
ponding  values  of  pressure  coefficient  will  likewise  be  different. 
Consequently,  Eq.  (93)  nay  be  separated  into 
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where  the  subscripts  A  and  i  refer  to  uie  upper  and  lower  surfaces. 

For  a  right  circular  cylinder,  Eqs.  (93)  and  (95)  simplify 
considerably  since  SN  and  r'-irfrx  are  zero.  Hence, 


where,  ar  usual,  ^  =  TTI^/4  and  c  =■  D 
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As  another  example,  consider  a  right  circular  cone  at  zero  angle 
of  attack.  For  tnis  case,  S/si  =S«*>st\S  a-id  for  reasonably  large 
'Sh  (say  3.0),  the  approximate  "reflected”  pressure  coefficient  is 
sufficiently  accurate.  Then  from  Eq.  (15)  or  Eqs.  (95), 
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which  when  substituted  in  Eq.  (9ir )'  yields 
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Kcw  that  the  zero  angle  of  attack  pitching  moments  have  been 
formulated  for  two  simple  satellite  shapes,  the  effect  of  large  tempera¬ 
ture  differences  on  the  equilibrium  angle  cf  attack  may  be  determined. 

For  the  case  at  hand  (sun- satellit e-Barth  system),  radiation  heat  tra  isler 
techniques  are  available  (Ref.  10)  for  approximating  the  upper  and  lc  ;r 
skin  temperatures.  For  example,  assume  the  following  average  temperature 
levels: 


T«,£  940°H 


Tw<  S  5I3°& 


(9?) 
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a: id  let 


T„  =  250 o°R. ") 
S~  =74  | 


(alt 


rfa) 


(100) 


rfhere  the  large  surface  temperature  difference  is  obviously  conservative* 
due  to  the  negligence  of  heat  transfer  within  the  vehicle.  Furthermore, 
for  stall  angles  of  attack,  the  static  pitching  moment  coefficient  is 
expressible  as 


which  for  equilibrium  conditions  (  Cvr\  -  may  be  re'-’ritten  as 


(101) 


The  pitching  moment  slope,  Cm*  ,  has  been  determined  for  the  cylinder 
and  cone  in  Sections  3. 2. 2.1  and  3. 2. 2. 3  respectively  for  the  case  of 
constant  surface  temperature.  These  equations,  however,  are  easily 
modified  to  include  the  ternperature  difference.  Consequently  for  the 
cylinder,  Eq.  (6U)  and  Eq.  (96)  combine  to  give 


(102) 
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(  radians) 


^Recently  published  results  (Ref.  11;  of  the  Discoverer  program  have  indi¬ 
cated  that  the  upper  and  lower  temperatures  are  quix-e  low  ana  not  much 
different  from  one  another.  In  some  cases,  these  temperatures  were  below 
freezing.  It  must  be  remembered,  however,  that  the  actual  equilibrium 
temperatures  achieved  are  strongly  dependent  upon  the  emissive  and  abs'rb- 
tive  properties  of  the  surface. 
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uhi re  the  values  of  Eos.  (9?)  and  (HOD)  have  been  used.  It  if  appaient 
iron  Eq.  (102)  tl.a'.  the  equilibrium  angle  is  very  small,  as  expet  tec,  for 
a  cylinder  with  a  large  temperature  difference  between  the  upper  .and 
lower  surfaces.  It  rust  be  remembered,  however,  that  the  temperature 
elfect  on  reflection  coefficients  has  not  seen  considered  in  this  analysis. 

The  cone  equilibrium  angle  may  now  be  found  an  a  similar  manner. 
Using  Eq.  (73)  and  Eq.  (98),  Eq.  (101)  becomes 
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where,  as  before,  Eqs.  (99)  and  (100)  are  used.  For  a  representative 
x**JJL  of  2/3  (foi  simplicity),  Eq.  (103)  is 
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which  is  very  small  for  reasonable  cone  angles. 

Although  this  analysis  is  obviously  incomplete,  it  may  be  qualita¬ 
tively  concluded  that  temperature  differences  between  opposite  sides  of  a 
satellite  have  little  effect  on  the  aerodynamic  pressure  difference,  and 
hence,  little  effect  on  the  equilibrium  angle  of  attack.  It  should  be 
remembered,  however,  tnat  although  solar  radiation  was  considered  herein, 
the  effect  of  solar  photon  pressure  on  satellite  pitching  moments  was  not 
included.  If  desired.  Appendix  A  may  be  used  to  calculate  the  additional 
moment  contribution  due  to  solar  pressure.  A  brief  check,  however,  has 
shewn  that  thi:  contribution  (for  a  cylinder  at  1  million  feet)  is 
smaller  than  the  aerodynamic  Cm0  due  *°  temperature  variations. 

3.3.L  Effect  of  Longitudinal  Temperature  Distribution  on  Stability 


It  may  be  recalled  from  previous  discussions  that  the  surface 
temperature  of  a  vehicle  in  free-molecule  flow  affects  the  surface 
pressure  (through  the  reflected  molecules),  and  hence,  aerodynamic 
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stability  or  the  vehicle .  It  v<is  shown  that  for  large,  conservative 
temperature  differences  te  tween  the  upper  and  lower  surfaces,  the  effect 
o-i  zerc-lift  pitching  moment  and  corresponding  equilibrium  pitch  angle  is 
very  small*  Another  temperature  effect  to  be  considered  is  the  influence 
of  a  longitudinal  temperature  distribution  on  the  pitch  characteristics 
of  a  satellite*  This  distribution  nay  be  the  result  of  (say)  conduction 
due  to  certain  ’’hot  spots"  within  the  vehicle,  free-molecule  convection, 
and/or  solar  radiation.  If  all  these  nodes  are  considered,  either 
separately  or  together,  the  resulting  temperature  distribution  becomes  a 
complicated  function  of  numerous  variables;  i.e.,  surface  material,  free- 
stream  conditions,  surface  orientation  angle  with  respect  to  the  sun,  etc. 
For  example,  if  a  radiation  heat  balance  is  considered  alone,  the  tempera¬ 
ture  ratio  T v/Te  depends  on  the  fourth  root  cf  a  function  of  the  angle 
between  a  surface  element  ard  the  sun.  The  square  root  cf  the  temperature 
ratio,  as  used  in  the  aerodynamic  pressure  equation,  then  depends  on  the 
eighth  root  of  the  sun  angle  function. 

The  laborious  details  of  a  heat  transfer  analysis  do  not  seem 
justified  for  the  rather  general  cases  considered  herein.  It  will,  there¬ 
fore,  be  assumed  that  an  average  temperature  distrioution  nay  be 
represented  by  a  simple,  axisymmetric,  power  law  expression  for  the 
purpose  of  weighing  the  effects  of  temperature  variations  on  the  aero¬ 
dynamic  pitching  moments.  For  example,  let 


where  ©  is  some  reference  temperature  ratio,  */£  is  the  non-dimensional 
body  length  measured  frer.  the  nose,  and  the  exponent  i)  is  a  constant. 

Now  for  a  body  of  revolution,  the  temperature  contribution  to  the  initial 
pitching  moment  slope,  Cm*  >  =37  be  found  from  Eq.  (25)  as 
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idler©  the  reference  area  and  length  areTE^/4  anrf  P  respectively. 
Now  for  a  right  circular  cone,  r  =  x-tom  ,  and  tr»  substitution  cf 
Eq.  (105)  into  Eq.  (106)  and  subsequent  integration  yield 


Since  realistic  v<*lues  of  ^  cannot  be  fai  different  from  zero  (say 
-14  ^  <  i  )j  Anri  the  two  se-ies  in  Eq.  (107)  tend  to  cancel  each 

other,  it  iaajr  be  seen  that  the  constant  temperature  (  ')  =  o)  increment 
oi  pitching  moment  slop which  is  small  when  compared  wiv.  tne  total 
patching  slope,  is  not  altered  appreciably  by  the  example  temperature 
distribution.  As  an  illustration,  Fig.  27  presents  a  ploy  oi  total 
Cm<  versus  the  exponent  ^  .  The  ^  s  o  conditions  are  noted  on 
the  figure. 

Other  example?  nay  be  performed  yielding  the  same  general  result, 
i.e.,  the  effect  of  temperature  variations  over  the  surface  of  a  satellite 
in  free-raoiecule  flow  do  not  significantly  affect  the  aerodynamic  forces 
and  moments  as  long  as  the  reflection  coefficients  tT  and  <rf  are 
assumed  invariant  with  temperature.  The  effect  of  variations  in  (T 
and  CT1  will  be  consicered  ii  the  next  section. 
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Figure  27.  Effect  of  Longitudinal  Temperature  Distribution  on 
Cone  Pitching  Mo"«:i\  Slope 
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3 . 3.5  Effect  of  Asy.Tn-etric  Reflection  Coefficient  Variations  on 

Stability 

As  Aontioneo  before,  the  effects  cf  variations  of  the  surface 
reflection  coefficient*  <r  and  <r  '  ,  are  analyzed  for  the  special  ca^e 
of  constant  surface  temperature;  in  particular,  for  Tw  ■  550°R.  For 
the  present  problem,  it  will  be  assumed  that  the  upper  and  lower  poi  tions 
of  a  body  of  revolution  have  different  values  of  the  reflection  coeffici¬ 
ents  C*  and  <T*  .  The  resulting  difference  between  the  upper  and  lever 
3hear  and  pressure  forces  will  then  produce  changes  in  the  pitching 
ncment.  The  magnitude  cf  this  asymmetric  loading  may  then  be  analyzed 
by  calculating  the  equilibrium  angle  attack,  which  is,  to  thr  first  order 
and  for  snail  angles,  defined  by  Eq.  (1'Jl/. 

Using  the  approximate  free  molecule  theory  (Ref.  h ) ,  it  may  be 
shewn  tnat  the  static  pitching  moment  coefficient  for  a  convex  body  of 
revolution  is 


which  is  the  sane  as  the  last  of  Eqs.  (25)  except  that  0“  and  (T*  are 
now  assumed  nependent  on  the  azimuth  angle  $  . 

Now  if  (T  and  (T ;  are  assigned  different  values  for  the  upper 
and  lower  surfaces.  Eq.  (10S)  becomes,  for  zero  angle  of  attack, 
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lir.ev'ifs,  the  ir  t’.al  pitching  nonent  slope  nay  be  found  by  differentiating 
Sq.  (106)  with  respect  to  erf  ,  »nd  taking  the  limit  as  o<  approa-hes 
zero.  The  result  is, 
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For  a  right  circular  cor.c,  r  =  Sir  - Tf £?/4  ^  c  =  C  ard 

Eqc.  1109)  and  \110)  integrate  to  yield, 
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where 


The  com  equal  vbi  jum  angle  of  att&cx  may  no,r  be  found  by  substitutirg 
Eqs.  (Ill)  sr.d  (112)  into  Eq.  (101);  i.e. , 
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For  purposes  of  representing  (at  extreme  case,  let  the  lower  surface 
reflect  in  a  ccmplete  specular  fashion,  0~i  =  CT*'  -  o  ,  while  the 
upper  surface  reflects  in  a  complete  diffuse  way,  (Tx  ■=.  =  1-0 

Then  A(T  -  AT'  -  -JL.  C  and  a  CT#1  -  '/z  •  Furthermore,  let 
S« -74-,  fTw/^tToT  =  O  (from  Fag.  2),  -ton£  “  </fs 

ana  xo**/£  st/j  .  fcltri  these  values,  Eq.  (113)  ^educes  to 
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The  value  of  Eq.  (11!*)  obviously  represents  an  improbable  condition. 
Otr.er  more  practical  reflection  coefficient  and  gecnetry  •'onbir.ations, 
however,  n^y  result  in  equilibrium  angles  of  attack  which  nay  signifi¬ 
cantly  affect  the  orientation  of  a  satellite.  Body  shapes  other  than 
the  cone  may.  of  course,  be  analyzed  for  the  effects  of  nen-uniform 
reflection  coefficients.  In  particular,  the  cylinder  is  very’  sensitive 
to  asymmetric  loadings,  ^specially  if  the  ends  produce  a  large  C**i0  • 

numerous  other  problems  nay  be  devised  which  involve  variations 
of  surface  temperature  and  reflection  coefficients.  Only  a  few  of  the 
no  re  extreme  (and  simple)  cases  have  been  discussed.  For  example,  it 
nss  been  assumed  that  <T  ana  0*'  art;  independent  of  surface  temperature, 
wmle  in  reality,  it  is  suspected  that  a  certain  dependence  between  the 
reflection  pher.oer.ena  and  surface  temperature  exists. 

3-3-6  Effect  of  Kon-Diffuse  Reflection  foeff i^ients  -  Stability 

In  this  section,  the  pitchirg  moment  slope  Cm*  is  investigated 
for  the  case  of  constant  surface  temperature  (550°R)  with  uniform,  but 
unequal,  values  of  the  tangential  and  normal  momentum  exchange  coeffici¬ 
ents,  o*  f.nd  <r<  . 
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Since  the  diffuse  -  '  )  values  of  have  beer, 

pi-esented  for  several  configurations  in  Section  3»2,  it  is  now  convenient 
tc  use  the  ratio,  Cm-t/Cmjtj,  ,  where,  in  general,  Ot.«  is  for 
c  4  (r‘4  i  ,  and  where  Cm^p  is  for  T  -  r‘  •=  i  .  The  examples 
chosen  for  illustration  are  the  cylinder,  sphere  and  cone. 

3. 3.6.1  Cylinder 


From  Eqs.  (62)  and  (6u)  it  is  easily  shown  that  the  cylinder 
0*y  ratio  is 
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which  is  mdependent^of  the  center  of  mass  location.  Now  if  the  approxi¬ 
mations  fTw  /  pfZ  =  O.Ob?  and  =?  7  4  are  used,  and  if  a 
representative  fineness  ratio  of  1.08  (an  optimum  value,  as  shown  later) 
is  assumed,  then  Eq.  (115)  nay  be  written  approximately  as 


„  ov?o  +  io8zr-o.c?i  <r' 

— -  —  — - — - —  vllo) 

CwyD  {  2S[ 

where  other  practical  choices  of  S«  ,  Tw  and  91 D  do  not  appreciably 
affect  tr.e  numerical  results.  Sq.  (116)  is  shown  plotted  in  Fig.  28 
versus  the  tangential  reflection  coefficient  (T  with  the  normal  reflec¬ 
tion  coefficient  <r‘  as  a  parameter.  The  ratio  C'/(T  could  also  be  used 
as  an  effective  parameter.  If  a  reasonable  lover  lisdt,  say  0.°,  is 
assigned  to  T  and  <T'  ,  it  may  be  seen  that  a  deviation  as  much  as  \8£ 
in  On*  could  be  incurred.  As  expected,  the  cylinder  Cm*  is  strongly 
aifected  by  the  shear  term,  and  hence,  decreases  substantially  with 
decreasing  \T  . 

3. 3.6.2  Sphere 


Section  2.1;  **■*  used  to  find  the  Cm*  io  fer  the  sphere 

since  Cm<  for  this  cor. figuration  is  directly  proportional  to  the  drag 
coefficient.  Then  from  Eqs.  (U?)  and  (51),  Cm*  may  be  written 

exactly  as 
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where  frcn  £q.  (50), 


Q{s.^piVe'$%  1W£-I  ^ 


Nov  if  the  usual  assumptions  of  So»  -  1  4  and  iTwA^frU  £  0-063 
are  used,  then 

r  2.  036  -*■  f-01 8 v” —  0 .94-<HT’ 
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which  is  graphically  presented  in  Fig.  29.  The  sphere  is  seen  to  be 
different  from  the  cylinder  in  that  o-  and  v  ’  nave  nearly  equal  and 
opposite  effects  on  Ct*.*  .  That  is.  decreasing  <r  ter.d3  to  decrease 
stability,  while  decreasing  cr1  tends  to  increase  stability  by  nearly 
the  sane  amount.  If  a  lower  iirit  of  0.8  i:>  assigned  to  the  reflection 
coefficients,  the  maximum  deviation  An  Cm**  is  again  seen  to  be  approxi¬ 
mately  18£.  For  <T=  cf'  •  however ,  the  deviation  is  very  snail.  In  fact, 
the  completely  specular  Cny^(T*T'  -  O  )  is  nearly  9$.?%  of  C*vp  . 

3. 3.6. j  Cone 

For  sen i-vertex  angles  greater  than  (say)  10°,  the  cone  pitching 
norent  si cue  at  zero  angle  of  attack  has  been  written  approximately  by 
So.  (73)  for  diffuse  reflection  and  by  Sq.  (11?)  for  general  values  of 
the  reflection  coefficients.  The  quotient  cf  Eq.  (73)  into  £q.  (11?)  is 


,\  .  (T  ff  ffW To*  1[.  1  xz^.  *A. 
Gn*  „|z(2-r-f  Jt  ip  2  xM 
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.  it- 


N-'V  ►  is  chosen  as  cex~'  fz  ,  then  =  */?  ,  and  the  center 

o:  mass  terirs  corr'enitnUy  cancel.  Eq.  (119)  tnen  reduces  to 

<±±  s  ±_±r-l-?£l£'  (120) 

where  r ;/  £*>  fU  as  again  assumed  O.O63.  Eq.  (120)  is  shown 
plotted  an  Fig.  3C.  Unlike  the  cylinder  and  the  sphere,  the  cone 
is  seen  to  be  affected  more  by  the  pressure  forces  than  py  tne  shear 
iorces  (at  least  for  the  low  fan* ness  ratio  core  selected  here).  An  18^ 
naxinur.  deviation  between  On*  at  <T-=l  O  f  (T*  -O  8  and  vice  versa,  as 
again  (by  coincidence)  evident.  The  cone  C**cx.  for  (J~ -  (T1  -  O-tJ  » 
however,  as  increased  by  only  6%.  From  Eq.  (120)  the  specular 
is  approximately  30^  higher  than  the  aiffuse  Ct*o^  * 


3-h.l  Review  and  Introduction 


In  Volume  I  of  this.  series  of  reports  it  is  emphasized  that 
successful  aerodynamic  stabilization  of  axisymretric,  rear -Earth  satel¬ 
lites  depends  upon  the  magnitude 
staoility  parameter  (pP*  .  This 


<J?  =  P 

■there 

y  -  Cmk  Stt  c 

I 

ft  *  atmospheric  density  at  c^-i-latus  rectum 
•>  »  perigee  radius 

Vft/V  ®  relative -to-inertial  velocity  ratio 

I  =  principal  moment  of  inertia  about  the  y  body-fixed  axis. 

It  is  also  noted  in  Volume  I  that  (?p  is  used  in  combination  with  the 
gravity  torque  parameter  ,  which  is  aefined  as 


(and  sign,  of  course)  cf  the  aerodynamic, 
parameter  is  defined  as 


vn 


(121) 


Ms  I  (123) 

1 

wnsrs  1 *  is  the  principal  nosent  of  inertia  about  the  X  body  axis.  In 
particular,  for  circular  orbits  it  is  shown  that  *ha  dominant  stability 
tern  (and  "spring  constant")  is  O*  +  2N\  . 

If  aerodynamic  stabilization  is  to  be  fessible,  tnen  the  question 
now  becomes  cr.s  of  determining  the  attainable  range  of  C 5p*“  for  any 
reala r tic  satellite  shape.  For  a  specified  orbit,  (?*  becomes  a  function 
of  1’  only,  and  hence,  the  problem  may  be  restated  as  being  one  of 
determining  th<;  ortiaua  proportions  of  a  given  general  shape  (cylinder, 
cone,  etc.)  which  make  the  stability  parameter,  P  ,  a  maximum. 


V 


For  the  configurations  considered  henceforth,  favorable  centor  of 
nass  locations  are  obtained  by  concentrating  a  homogeneous  payload  in  the 
forward  portion  of  the  satellite.  The  aft  oortiorsof  the  models  are 
assumed  to  be  rigid,  but  masiless,  shells  or  skirts.  It  is  obvious  that 
these  assumptions  ar-  not  realistic,  but  >t  i,.  hoped  tnat  the  design 
anal/ses  to  follow  v  11  serve  as  guidelines  for  the  more  sophisticated 
ces. gns  oi  aeredynam rally  stabilizea  satellites. 

Although  many  optimization  procedures  are  availaole,  the  technique 
used  herein  is  to  differentiate  T  while  holding  payload  mass  (or  weight) 
and  center  of  mass  cons  van  t.  Nov;  if  the  payload  density  is  specified, 
then  the  payload  volume  is  also  known.  Consequently  is  convenient  to 
non-dimensionalize  Eq.  (122)  by  multiplying  both  sides  oy  the  payload 
mass,  nv  ,  while  dividing  both  sides  by  the  payload  volume  to  the  one- 
third  power,  v ;  i.e. . 


pi  -C»wT7D*  m. 

/vy*  “  4Xv'/j 


(12b) 


where  tne  right  hand  side  is  independen?  of  m.  since  I  is  a  linear 
function  of  mass.  Accc.dingly,  Eq.  (12b)  is  a  function  of  satellite 
external  geometry  and  the  center  of  mass  location.  For  purposes  of 
comparing  several  configurations,  however,  the  center  of  mass  location 
may  be  conveniently  written  in  terms  of  volume  ratio;  i.e.,  the  ratio  of 
payload  volume  to  total  satellite  volume,  including  the  skirt. 

Since  optimum  geometry  is  to  be  deterrdrea  from  the  maximization 
of  r(m/v^)  ana  since  the  center  of  mass  is  specified,  the  gravitational 
parameter  M  is  calculable  and  is  presented  m  the  following  analyses  of 
the  various  configurations.  Also,  it  should  be  noted  that  the  maximiza¬ 
tion  of  T  tends  to  minimize  X  ,  thereby  tending  to  maximize  M  • 

Tne  oody  chapes  chosen  for  the  following  design  studies  are  the 
same  as  treated  in  Section  3.2.  The  basic  body  shapes  -  cylinder,  sphere 
and  cone  -  are  analyze a  separately  and  critically  compared.  The  cximun 
r's  of  the  power  law  body  and  the  truncated  cone  are  snown  to  be  not 
much  different  frea  those  for  the  pointed  cone-  and  consequently,  these 
bodies  are  not  extensively  compared.  The  hemisp.  re-cylinder  ard  r*one- 
cyiinder  are  compared  with  the  basic  bodies  vhicn,  when  taken  together, 
ma.e  these  composite  configurations. 
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Figure  31.  Cylinder  Nomenclature 


’Iflglecting  the  contribution  of  the  skirt  of  the  vehicle  to  its 
.non/Iflt  of  inertia  about  the  y  axis,  I  ie  found  to  be 


(125) 


and  for  later  use,  the  mcsn>»nt  of  inertia  about  the  X  axis  is 

r  «  s®*  •  ass) 

8 

For  a  cylindrical  homogeneous  payload 

=■  t/z  (127) 


vhere,  of  c«*rse. 


v  =  Hlfrl/V 
4  l'?' 


^■y*  .  fVv/  \  TT  in5* 

D  '  [1  '"M  4.  DJ  '  (128) 

S"Srs,«?“'"g”K“s,;"’  *■ -  v  ™ 

=  74rf°i  fcr, '  “4  i'-.'vVl 

~ - - - - - i  .  (12?) 

SST**”*  Eq-  (12?)  *«*  to  i/o  1  setting  equax  10  «ro 

s^-sriS  IrlHSSKI 

becomes  approximately,  J  ±'J'  t/*‘  *  therefore 

‘*4L-‘RS8L-^.o  <» 

Which,  for  volume  ratios  of  1/1, ,  1/2  j/i  ,ntf  ,  . 

solutions:  /  >  /  »  ->/*♦  ana  i,  has  ~he  graphical 


V/VV=!/t  :  lWoU=4Si 

^  -k  :  (Vd\w  s  /80 


V/VT  =  V*  ;  sAi o)eM.=  1.08 

<K  =  i  :  (Wo)  s  0  7fo  • 


Siibstituting  the  foregoing  results  into  Eq.  (129)  gives  the  sawttun  values 
ol  .  Moreover,  .f  Uie  payload  density  is  specified,  r.** 

nay  be  round  as  i  functior.  of  payload  mass  or  weight.  The  following  table 
sumarizes  these  results  for  the  case  of  Vx.  slua/ff  a  1_  AJS0 

sheen  arc  the  values  of  the  gravitational  parameter,  .-.i.  -  /T,/l'\  .a 
which  correspond  to  the  optimum  fineness  ratios.  ”  has  tne  dimnsicns 
cf  ft/slug.  Also,  the  payload  weight  W  is  evaluated  at  sea  level  for  a 
nean  acceleration  of  gravity  of  3?.?  ft/sec11. 


Table  1.  Cylinder  Sunnary 


[y'0)oer  l(r?v.) 


ion 


The  resui.s  shown  in  Table  1  (except  for  Y/vy  =i  o  )  are  presented  as 
plots  of  r  versus  payload  veight  in  Figs.  35,  36,  and  37.  The  signifi¬ 
cant  features  of  these  results  art*  <*«  Section  3. I*. 2  2* 


^The  payload  de  isity  of  1/2  sltg/ft3  represents  a  reasonable  average  of 
U.S.  sate1 lite  payloaos  launched  to  date. 
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v/vT  =  ^(s-'Ze'i 

-fy>  -  ,  (135) 


Hcreover,  the  moment  of  inertia  about  the  X  axis  is. 


T  -  S  1 

x  " 


which  vhen  combined  with  Eq.  (133)  leads  to 


j  +2€T“4r^'> 
Tit  +  #<=*• -JL? 

C?  ?  IO 


(136) 
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Substituting  Eq.  (69)  for  a net  Eq.  (133)  for  I  into  £q.  (12U)  yields 


r  - 

i  — —  j 

\aCc  \  .  -2Jl  ’ 
l  T  5S*. 

sw-vytf 

fr\v5  £  | 

V*r)  (5-Z6)1’ 

(138) 


where  the  geonctrical  terms  are  all  functions  of  or/VV  since  x<*«./D  = 
-fC<0  -  -f  (ir/VT\  .  The  term  in  the  brackets  may  be  estimated  using 
Eq.  (50)  and  Fig.  2,  and  consequently  has  a  value  near  2. Hi.  Although 
Eq.  (138)  may  be  written  in  teres  of  -v/Vy  ,  the  lengthy  expression  is 
quite  cumbersome . 

It  .-ore  frT  .ble  to  construct  plots  of  x«*»/D  and  'v/'Vt 
versus  the  depth  ratio  ^  .  Accordingly,  these  curves  are  shown  in 
Fig.  33  along  with  variation  of  the  gravitational  parameter  M  .  It  is 
noted  that  Cor  all  possible  values  of  fc  ,  M  ^.0  and  hence  the  gravita¬ 
tional  •torques  always  accent  the  aerodynamic  torques  in  the  case  of  a 
"loaded*  sphsie. 
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Table  2  presents  a  summary  ox  the  aerodynamic  parameters 
developed  above.  Here  again.  the  payload  denr ixy  is  assumed  to  be 
i/2  slug/rt3.  Discussion  cf  the  results  ana  comparisons  with  the  cylinder 
and  cor.e  are  noted  in  Section  3-U.2.U  and  Fig3.  35,  36  and  3?. 


Table  2.  Sphere  Stewary 
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Figure  SU.  Cone  Nomenclature 
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Furthemcxt  it  is  easy  to  show  *hat 


which  leads  to 


IA  =  tliffi  . 

4-  +  «H-zc, 


Oil) 


Following  the  procedure  used  for  the  cylinder  and  sphere,  Eqs-  (73), 
(139),  and  (UiO)  nay  be  substituted  into  Eq.  (121i);  i.e., 


The  question  now  becomes  one  of  determining  the  optimum  cone  angle  S 
which  rendeis  Eq.  (1U2)  maximum  for  a  given  volume  ratio.  The  temperature- 
dependent  term  in  Eq.  (11:2)  obviously  complicates  the  maximization  process. 
As  nay  \n  seen  in  Fig.  25,  however,  this  term  is  snail  when  compared  with 
the  first  term,  and  consequently,  contributes  very  little  in  the  ue ter¬ 
mination  of  an  cptimua  cone  angle.  Then,  by  neglecting  the  temperature 
term  and  differentiating  Eq.  (2:2)  with  respect  to  £  ,  and  setting  equal 
to  zero  yields, 


Z$U=  ^Vr  =  ,Z 


iifi3) 


which  is  independent  of  volume  ratio.  The  substitution  of  Eq.  (11:3)  into 
Eq.  (Hi9)  gives  the  maximum  (nearly)  value  of  Pm/fv»  5  !•«•» 


(Vm/v**)  = 

'  ''fcvx'r 


Using  Fig.  2  for  a  representative  value  of  )[Tw/T»  ,  Eq.  (lijli) 

reduces  to 

(r^L  =  ?5.s  [i - ?4(^} /(^t) .  {nW) 

Furthermore,  the  optimum  cone  angle,  Eq.  (Hi3),  nay  b«  s"*bs  Mtuted  in 
Eq.  (liil)  to  give 

Nl  =  '/■}  (1-6) 

*-hich  is  positive  and  stabilising.  The  foregoing  results  are  summarized 
in  Table  3  and  Figs.  35  through  33  for  volume  ratio?  of  l/U,  1/2,  3/1 
and  1  and  for  the  case  of  m/-v  -  i'z.  siujiVf* 


Table  3.  Cone  Summary 
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1  iA 

35.3° 

O.U73 

75.5 

95.2 

96U 

0.333 

i 
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36. 

369 

0.333 

§ 

35.3° 

0.681 

15.75 

19.2 

1911 

0.333 

1 

->5.3° 

0.750 

8.95 

11.3 

HU 

0.333 

3-i.2.!i  Comparison  of  Basic  Body  Shapes 


Now  that  the  stability  parameters  of  three  basic  satellite  shapes 
have  been  calculated,  it  is  desirable  to  analyte  the  advantages  and  dis¬ 
advantages  of  each  configuration  from  the  point -of -view  of  maximizing 
stability  for  a  specified  payload  condition.  The  principle  used  in  the 
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preceding  sections  is  to  optimize  I1  while  holding  the  payload  volume 
(or  mass)  constant,  and  then  observe  the  variation  of  this  maximum  as  the 
center  cf  mass  and  payload  veight  of  each  configuration  are  varied.  In 
otner  words,  the  question  may  be  asked,  "Which  configuration  has  the  most 
aerodynamic  stability  for  a  given  payload  veight  and  center  of  mass?" 

This  form  of  optimization,  however,  does  not  completely  solve  the  problem 
in  that  maximum  |*  for  a  specified  payloai  volume  (cr  mass)  and  center  of 
mass  also  specifies  tne  diameter  of  the  satellite.  Consequently,  the 
preferred  configuration  from  a  stability  viewpoint  may  have  a  prohibi¬ 
tively  large  diameter. 

The  results  of  the  stability  analyses  of  the  three  basic  body 
shapes  are  shown  plotted  in  Figs.  35  through  38  for  the  four  volume 
ratios.  Since  fully  loaded  homogeneous  cylinders  and  sphere  are  neutrally 
stable,  only  the  cone  appears  in  Fig.  38  for  v/'Vr  -  l  O  -  Also  shown  are 
integer  values  of  satellite  diameter  along  each  maximum  f  line.  In 
Figs.  36  and  37,  values  of  V  for  constant  diameter  are  shown  to  illus¬ 
trate  the  deviation  from  the  optimum  T  .  These  curves  clearly  show  the 
difference  between  maximum  T  for  constant  diameter  and  maximum  P  for 
constant  mass  (except  for  sphere).  It  may  be  seen  from  Figs.  35,  36 
and  37  that  the  cone  has  the  most  aerodynamic  stability  for  a  given  pay- 
load  weight,  while  the  sphere  has  the  least  or  nearly  the  least.  The 
cone,  however,  also  has  the  largest  diameter  for  a  given  payload  weight. 
For  example,  consider  V/vY  =  '/z.  and  a  weight  of  300  lbs.  (Fig.  36). 

The  optimum  cone  is  about  6  feet  in  diameter  while  the  cylinder  is 
only  3  ft.  The  sphere,  which  has  the  least  suability,  requires  a  h  ft. 
diameter. 


Now  suppose  there  is  a  requirement  that  the  satellite  be  no 
larger  than  3  ft.  in  diameter.  Then  at  these  conditions  (W  -  300  lbs  , 

■v* |vT*Vi  f  D  =  34*.)  the  cylinder  now  appears  the  most  stable.  In  fact, 
the  2  ft.  cylinder  is  better  than  the  3  ft.  cone.  At  Y/Vy  -  V4  and 
V/  =  300  lbs.  (Fig.  37)  the  optimum  cone  has  a  diameter  near  5  ft.,  while 
the  cylinder  is  still  neax  3  ft.  The  sphere,  which  is  slightly  more 
stable  than  the  cylinder  for  this  volume  ratio,  is  about  3-1/2  ft.  in 
diameter.  If  a  3  ft.  minimum  diameter  must  be  held  for  the  v/'Vy  -  3/4- 
case,  however,  the  cone  has  the  irost  stability.  If  the  length  of  the 
satellite  is  no  problem,  this  example  may  be  carried  still  further.  That 
is,  for  Y'/v^.  -  i/4.  (Fig.  35)  the  300  lb.  payload  may  be  carried  by  a 
cylindrical  satellite  of  less  than  3  ft  diameter.  Furthermore,  the  value 
of  P  for  this  case  is  considerably  higher,  as  expected,  than  for  the 
other  cases  with  higher  volume  ratios. 

A  few  words  are  appropiiste  a'oout  the  combineu  stability  due  to 
aerodynamic  and  gravitational  effects;  more  specifically,  the  quantities 
OV*  and  3M  .  Because  of  the  basic  difference  between  (?p*  and  3M  , 
a  combination  cf  these  parameters  is  impractical  to  optimize  analytically. 
O*  ,  he. .ever,  is  a  function  of  T  and  the  orbit  conditions,  and  there¬ 
fore,  maximum  P  maximizes  (?p*  for  a  given  orbit.  Furthermore,  it  is 
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desirable  to  keep  M  positive  (low  fineness  ratio)  to  auvuent  the  aero- 
dynamo  stability  in  the  examples  above,  or  la  Pigs.  35736  and  37, 

SDher,e  ®!d  »■*«  positive  values  of  f-1  ,  whale  for 

“n  op. ui.ua  cylinoer,  M  _s  near  zero  or  negative.  Consequently,  before 
e?eeotleCe=  3  satei‘-lte  Shape  based  on  naxir-un  aerodynanic  stability,  the 
ef.ects  of  gravitational  torques  should  be  considered.  As  shorn  before 
o^a^?^e'X'CeR  ae,r?^®ic  “<*  gravitational  effects  depends  on  the 
orbit  tonomons  as  veil  as  the  satellite  shape,  and  that  naxuaun  aero- 
djmar.ic  stability  certainly  does  not  necessarily  prevent  unstable  notions 
at  all  near-Earth  orbit  conditions.  These  latter  observationTl-e 
particular^  true  at  higher  altitudes  where  the  values  of  (?*  cry  become 
snail,  oving  to  the  decreased  atnosphenc  density.  ^ 
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Sc/9  Le*£L  Yf E.IQHT  •*  LffS 
Figure  33.  Ccesoarison  of  Basic  Body  Shapes  -  -  V+ 


5£/»  Levee  y/e/QHr^  Lbs. 

Figure  36.  Cenrarison  of  Basic  Body  Shapes  -  'V/Vj  ■*  Vx 
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3.11.3  Con?? ex  Body  Shapes 
3.ii*3*l  rower  Law  Body 

Following  the  analyses  of  previous  body  shapes,  the  static 
staoility  parameter  Pm/v’/*  as  well  as  V  and  M  will  be  examined  for 
the  power  law  body.  As  before,  the  body  will  be  assumed  composed  of  a 
nollow,  massless  skirt  with  a  homogeneous  payload  as  shown  in  Fig.  39. 


Figure  39.  Power  Law  Body  Nomenclature 


Tne  moments  of  inertia  of  a  power  law  body  may  be  determined  as 


I  -  rn&r?r'{7.n+‘) 

8  (4n+i)  (U*7) 

M  s  Lil  = 

I  4  (4mi)(!!d)z 


ana  tr.e  pertinent  geometric  relationships  are 


«<cvn-  _  aJ 

jL  "  2n*z 


V' 


yt  =  iro  ^ 


4<2«fO 


Yv*  W-\JH _ H\  i  3 

—  *  |UrTj4{>..Hl)  l  t>  )j 


'Jamg  the  approximate  pitohmg  none, it  slope,  Lq.  ( ?6),  in  which  the 
temperature  depeneort  tern  is  neglected,  the  ae-odynardo  stabilitv 
paras* tor  f  m  y-rVj  Jn  terns  0e  volume  and  mass  parameters  is 


r  = 

■fcffpwQ-  »V  [z' n  H)  ~  (MET  ImMH 

[(%p(2«sX«*)1  +  4(wi)ftf  ] 


The  optimum  fineness  ratio  for  maximum  stabilitv  for  a  piven  payload 
volume  or  miss  nay  be  found  by  differentiating  Eq.  (lhy)  with  respect  to 
XIC  ana  netting  equal  to  aero.  This  operation  yields 


(M 


and  ~her.  this  equation  is  substituted  into  the  last  of  Eqs.  (Ih7),  the 
inertia  parameter  M  beccnes 


which  is  independent  of  body  curvature  and  7olume  ratio. 

Substituting  Eq.  (150)  iiito  Eq.  (lh 9)  gives  the  approximate 
maximum  values  of  Pm/v^s  for  gj  yen  payload  volume  or  mass.  Plots  of 
Eq.  (Hi9),  using  Eqs.  (H48),  are  shown  in  Fig.  1*0.  It  may  be  seen  that 
the  value  of  the  exponent  has  little  effect  on  maximum  stability  for  gi^en 
volume  ratio.  It  should  be  remembered,  however,  that  although  the  curves 
in  Fig.  IjO  nay  oe  extrapolated  to  rv  =  o  (cylinder),  they  become  invalid 
at  “t  =  o  ;  i.e.,  Eq.  (liiO)  applies  only  for  o<ft£  i  .  The  exact 
cylinder  values  for  (  P m/Yv»)  ,**•<  are  shown  on  ilic  ordinate.  The 

difference  between  the  extrapolated  values  and  exact  cylinder  values  is 
due  to  a  discontinuous  shift  of  force  center  as  a  decreasing  (n-*o) 
power  law  boqy  suddenly  becomes  a  flat-faced  cylinder.  Also  shewn  are 
the  more  accurate  values  of  (rm/yV/m(y(  for  the  rt=l(conc)  and 
n  ~  Vi  bodies;  i.e.,  including  the  temperature  term  in  the  solution  of 

*a.  (75). 

It  also  appears  from  Fig.  hO  that  the  value  of  n.  whicn  gives 
maximum  pm /tTv*  is  something  greater  than  unity.  Since  n  greater  than 
unity  represents  a  body  with  a  concave  surface,  these  configurations 
cannot  logically  bp  treated  using  conventic-al  free -molecule  theory.  That 
is,  the  basic  assumption  of  free-molecule  theory,  as  apuiied  to  gasdynamic 
calculations,  is  that  the  incident  momentum  flux  is  not  disturbed  by  the 
presence  of  the  body.  Regardless  of  this  difficulty,  Vvever,  the 
"flatness"  of  the  curves  in  Fig.  I.C  indicates  that  the  *  itimum  body 
probably  does  not  have  much  more  stability  than  the  cone  ( n  =  l)  . 

Table  U  shows  the  various  maximum  stability  parameters  along 
with  (£/6)9pt  ,  and  Xor*H  as  calculated  from  Eqs.  (lu£),  (lh9)  and 
(150).  W  is  the  sea  level  weight,  m7‘>r  =  ,  and  P  is  in 

ft/slug.  From  Eq.  (151)  the  inertia  pa-aaeter  M  is  approximately  1/3 
for  all  power  lav  bodies. 


It  nay  be  noticed  ♦Mt  the  for  n  =  1  (cor.e)  m  Table  h. 

are  slightly  ( ~  U% )  lesa  than  the  cor.e  values  shown  in  Section  3*^.2. 3- 
This  nifference  is  again  due  to  the  negligence  of  the  tecperature  depen¬ 
dent  term  in  Eq.  (75)  for  Crv^  . 


Figure  ill.  Truncated  Cone  Nomenclature 


For  the  case  of  a  "loaded"  truncated  cone,  the  pitch  and  roll 
aonents  of  inertia  about  the  center  of  mass  are  functions  of  the  diaaoter 
ratio,  £  .  as  well  as  the  payload  length  ratio,  fl\  .  Consequently  it  is 
not  difficult  to  show  that 


I=  80  j  ) 
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x“  40 
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cTbImo? 


052) 


where 
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and  where 


w  1  ^ 


Likewise,  the  other  geometrical  properties  may  be  expressed  as 

r£  _  |  "fO'PJ  +3.-^(Vl1~fe)  fr  V  ' 

**T "  1  L  '  + 
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The  first  of  Sqs.  (1$6)  may  be  solved  for  ^  m  terms  of  b  and  the 
volume  ratio}  i.e.,  v 


(15?) 
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Eqs-  (15b)  md  (157)  way  be  used,  if  desired,  for  rewriting  the 
equations  for  ar.d  X  in  terma  o£  £  ,  *!&  ,  and  'r'/'Vr  .  The 

resulting  expressions,  however,  are  unwieldy  for  rapid  nunerical  calcu¬ 
lations.  Therefore,  the  stability  parameter  P  ,  as  defined  before,  will 
now  be  written  using  only  Eqs.  (80),  (1 5?)  and  (156).  The  resulting 
combination  yields 


For  given  volume  and  diameter  ratios,  the  approximate  optimum  fineness 
ratio,  £(D  ,  which  gives  maximum  stability,  may  be  found  by  differenti¬ 
ating  Eq.  (1§8)  with  respect  to  Xl 0  .  setting  equal  to  zero,  and 

solving  for  *(£)  .  The  resulting  egression  is 


m. 


~  !  h  c 


(159) 


Furthermore,  Eq.  U5?)  combines  with  the  last  of  Eqs.  (152)  to  give  the 
optimum  value  of  the  inertia  parameter  j  i.e., 

M  =  ’/?  (i6o) 

which  is  the  same  as  for  pointed  cones  and  power  law  bodies.  The  maximum 
value  of  r 3  may  now  be  found  by  substituting  Sq.  (159)  into 
Eq.  (156).  Ac  mentioned  alive,  C  ,  &  ,  Xc-r.i JL  ,  and  ^  are  only 
functions  ox  diameter  ratio  £  and  volume  ratio  'T/Vy  j  and  hence,  the 
maximum  Pm/VVb  depends  only  on  these  ratios. 

Fig.  1x2  presents  plots  of  Eq.  (158)  using  Eq.  (159)  for  various 
values  of  the  volume  ratio.  As  in  the  case  of  power-law  bodies,  the 
curves  arc  not  extended  to  the  cylinder  abscissa  (fc  =  I  O)  since  the 
equations  lor  Cm*  used  herein  are  invalid  for  cyl'  -ders.  The  exact 
cylinder  values  for  (P  however,  are  show,  on  the  ordinate. 

Moreover,  representative  values  of  the  more  exact  (P»w/v‘v*)m*x  are  also 
shown  using  from  Eq.  (79)  with  ^*=7  and  ( i/s~)  frwT«*  £  .065 

frem  Fig.  2.  It  may  be  seen  that  the  approximation  of  Eq.  (80)  is  in  very 
little  error.  It  also  appears  from  Fig.  Jj2  that  maximum  I"1 
occuis  for  configurations  with  non-zero  diameter  ratios,  which  is,  in  -act, 
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5  £  -gnificant  and  useful  development.  Then  for  a  specified  payload  volune 
cr  mss  requiremc.it,  the  pointed,  conical  satellite  ioes  not  represent  the 
mist  stable  configuration  as  originally  suspected.  The  blunted  cone  also 
represents  a  more  useful  snape  iron  the  viewpoint  of  re-entry-,  payload 
(.aikagiug,  and/or  antenna  uesign. 

The  results  of  the  stability  analysis  of  the  truncated  cone  are 
summarized  in  Table  5*  Here  again,  the  payload  density  is  taken  as  l/2 

sing/rO. 


Table  $•  Truncated  Cone  Summary 


loo 


It  should  be  re-empha sized  that  the  pointed  cone  values 
are  slightly  less  than  the  cone  values  in  Table  5  due  to  the  anission  of 
the  temperature  term  in  Sq.  (’>9)  for  CmQt  . 

Before  concluding  this  section  it  should  be  noted  that  there  is 
a  ri*d2?**<*y  between  the  results  of  the  power  law  bod/  and  the 

truncated  ccr.c.  In  particular  if  3  is  replaced  by  i-n  (  n  equals 
exponent),  the  first  two  terra  of  Taylor  expansions  of  the  body  profiles 
about  the  point  x/j l  =  1  are  identical.  This  similarity  may  be  noted  by 
corparing  the  values  of  ( f  m'r'/j)  m  for  a  given  'V/vv  in  Pigs.  J»G 
and  h2. 
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3-U-3.3  Hemisphere-Cylinder 


Figure  U3.  Hemisphere-Cylinder  Nomenclature 


As  in  the  problem  of  determining  the  pitching  moment  slope,  the 
calculation  of  the  inertia  end  geometrical  characteristics  of  the  hemi¬ 
sphere-cylinder  is  considerably  more  complicated  than  for  the  basic 
shapes.  The  parameters  shown  below  are  written  for  the  two  conditions  of 
tnc  payload  fineness  ratio  •€.  greater  than  and  less  than  1/2.  Since  the 
skirt  is  assumed  massless,  the  sphere  equations  of  Section  3.h.2.2  may  be 
utilized  for  £  <  Vz.  .  Consequently,  it  is  not  difficult  to  snow  that 
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(162) 


vhich  nay  be  conbined  to  give 

_  j  --§-£-+ 2 


M 


n  = 
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The  results  of  tn.*»  sphere  analysis  lead  to 


Xc«y  _ 
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Tne  total  fineness  ratio  and  payload  voluae  nay  also  be  written  in  terms 
of  the  payload  fineiess  ratio;  i.e.. 
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(165) 


(366) 


The  stability  parameter  P  and  its  alternate  forms  may  now  be  formulated 
using  Eq.  (fib)  for  Cvw*  and  the  foregoing  Eqs.  (l6l),  (16b)  and  (166). 

If  desired  Eqs.  (165)  may  be  solved  for  €  _  in  terms  of  Ifp  and  'y’I'Vj 
and  then  substituted  in  the  equations  fo’’  I,  x&m/D  and  Y" .  As  in 
the  cases  of  the  sphere  and  truncated  cone,  however,  the  resulting  expres¬ 
sions  do  not  necessarily  simplify  the  problem.  Consequently, 
is  written  here  as, 


1  m(X' 
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where  X  and  ^  are  defined  by  Eqs.  (85).  If  and  (  Vsu#)  fTWT*> 
are  representea  cy  the  typical  values  in  Fie.  2,  then  X  and  Y  are 
approximately 

X  =  o.WoWof  J  0 
Y  5  0  4«o(&>)  +  1-195 

where  it  is  understood  that  %tP  is  always  greater  than  1/2.  Moreover, 
it  should  be  re-emphasized  that  €  ,  in  Eqs.  (167),  is  a  function  of 
£<’0  and  *vY*A  -  For  references  the  solutions  if  Eqs.  (165)  ire 


Because  of  the  obvious  difficulties  involved  in  differentiating 
Eqs.  (167)  with  respect  to  ^Iq  ,  setting  equal  to  Luo,  and  solving  for 
(£/D)opt  >  n0  attempt  has  been  made  to  analytically  determine  tne 
maximum  value  of  pm/V^  for  henisphere-cj  binders  with  total  fineness 
ratios  greater  than  1/2.  The  stability  parameter  P  ,  however,  may  be 
plotted  against  payload  weight  for  constant  values  of  r n/v,  'WVT  and 
either  i(t>  or  0  .  Constant  X(P  curves  are  presented  in  Figs.  Wj, 

U5  and  h6  for  the  usual  i^/tt  -  <1  and  V/Yr  equal  1 /h>  1/2 

and  3/U.  The  examples  chosen  are  the  limiting  case  of  jf / 0  ~ which 
is  a  hemisphere,  and  XlO  -  Wc>  which  is  equivalent  to  the  ronplete 
sphere  with  the  same  payload  volume  ratio.  Curves  of  constant  diameter 
(=  2  f wi)  are  also  shown  in  Figs,  uii,  h$,  and  U6  for  the  three  volume 
ratios.  Since  varies  along  the  constant  diameter  curves,  a 

graphical  solution  for  maximum  Pw/vW  ,  and  hence  (X/0)o9r  ,  may  be 

possible.  Such  a  solution  is  shown  in  Fig.  Uh  for  in  which 

the  optimum  fineness  ratio  appears  to  be  near  I*. 5.  The  corresponding 
optimums  for  the  other  volume  ratios,  1/2  and  3/h,  are  not  apparent. 

Included  in  the  comparisons  of  the  hemisphere-cylinders  are  \  e 
basic  sphere  and  optimum  flat-faced  cylinder.  Of  the  configurations 
shoe*.,  it  is  clearly  seen  that  the  limiting  case  c.  ‘he  hemisphere 
(lfD -We)  is  the  most  stable.  This  configuration,  nowever,  is  probably 
impractical  since  the  diameter  requirement,  as  noted  in  the  figures,  is 
the  largest.  The  optimum  hemisphere-cylinder  (Fig.  fjl*  cnly)  is  seen  to 
be  only  slightly  better  than  the  optimum  flat-faced  cylinder.  It  is  also 


noted  that  the  equivalent  hemisphere-cylinder  (V/d7  sane  a?  for  sphere) 
has  a  riucn  higher  value  of  P  than  the  corresponding  sphere.  This  is 
er'pected,  however,  since  the  back  surface  of  the  sphere  contributes 
essentially  nothing  to  the  stability.  Furthermore,  it  is  evident  that 
the  sphere  with  r/Yf»(/4  (Fig.  lib)  and  the  herd. sphere  with  'W'V'r  '  */2- 
(rig.  IS)  have  the  sane  stability  (  e  is  the  sane  for  both  configura¬ 
tions.) 


tifym  Ji'i.  octpsriocn  of  Sphere-Cylinder  Configurations  -  v/v t  *  Va 
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no 


iJi.l.lx  Cone-Cylinder 

The  ccnc-cylinder  analysis  will  be  limited  to  the  case  where  the 
cons  is  completely  filled  with  payload  and  the  cylinder  is  empty.  The 
geometry  of  this  configuration  is  shown  in  Fig.  Ii7. 


Figure  U7.  Cone-Cylinder  Nomenclature 


The  principal  mOiaent  of  inertia  about  the  y  axis  is  easily 
found  using  Bq.  (139)  of  the  cone  analysis:  5*e., 


;=  wdT-,  *  JLfSaffXTl 
’  'go'  T^-f)  Vo)  j 


and  likewise  I*  is  evaluated  as 


X  =  JL-yvtlT  . 

4-c 

Sqs.  (3?C)  and  (171)  combine  to  give 
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For  the  hcn:geneous  conical  paylcad  sketched  in  Fig.  li?, 


9^) 

1  4-t  &(v/vr) 


(173) 


Using  the  cone-cylinder  Cm*<  from  Eq.  (88)  and  the  equations  above, 
Eq.  (12h)  for  f T*n/vv^  “7  be  expressed  as 


'-here  J  and  K  are  defined  by  Eqs.  {£?).  The  first  of  Sqs.  (173),  of 
course,  may  be  utilized  to  sake  Eq.  (17h)  independent  of  Xom /JL- 

As  in  the  case  of  the  hemisphere- cylinder,  the  maximization  of 
Eq.  (17b)  by  analytical  methods  is  not  practical,  it  .  possible,  .lcvever, 
to  plot  r  versus  payload  weight  (say)  for  constant  values  of  m/v^vY'vy 
and  D  •  The  optimum  fineness  ratio  may  be  feund  by  drawing  a  constant 
i/0  line  tangent  to  the  constant  diameter  curve.  The  optimum  x/D  line 
then  represci"?  the  cone-cylinder  configuration  which  has  the  maximum  P 
for  a  specified  payload  weight  and  volume  ratio.  This  graphical  optimiza¬ 
tion  procedure  is  illustrated  in  Figs.  U8,  h?  and  $0  for  a  payload  density 


of  1/2  slog/ft3;  volute  ratios  of  l/lt,  1/2  and  3/1;  and  a  diaraster  of 
2* fast.  Also  plotted  in  these  figures  are  tie  equivalent  curves  -or  -he 
opti-j*  cone  and  the  ootfauus  cylinder.  It  should  be  remembered  that  the 
cone-cylinder  configuration  used  in  these  comparisons  is  composed  of  a 
‘iilly  ioadea  cone  and  a  massless,  cylindrical  skirt,  ttien  ‘y^T-  /♦, 
it  is  seen  from  Fir.  US  that  the  optima*  cone-cylinuer  has  tl.e  lowest 
value  of  f  for  a  given  weight;  heweve. .  when  TfW  -  A  or  /4-  the^ 
optimum  cone-cylinder  has  a  higher  value  of  "  than  the  op-iBu*  . 

It  is  also  seen  iron  there  figures  that  for  a  given  diameter  and  volume 
ratt"  ther-  -a  a  range  of  sea-level  weights  for  which  the  cone--ylin-e. 
has' the  highest  value  of  P  of  the  three  configurations  shown,  ror 
v/Y-  w  l  O  .the  cone  and  cone-cylinder  are  obviously  coincident. 


Sea  Level  Weight  -  ‘t»s 

/iyyre  L9.  Cotepsriscn  of  Cone-Cylinder  Confirureticns  -  V/vT  =  */z 
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h.  PiKAKIC  STA8ILI1Y  ANALYSIS 


h.l  General  Considerations 


In  Volume  X  of  this  senes  of  reports,  it  is  argued  that  the  aero- 
(^T.rdc  damping- -in-pitch  of  near -Earth  satellites  produces,  for  all 
practical  purposes,  a  negligible  effect  on  the  satellite's  oscillatory 
notion.  Even  if  the  dating -in-pitch  derivative  Cm<y  is  considerably 
larger  than  its  static  counterpart  Cm*  »  the  non-dimensional  pitch  rate 
®£  /xV»,  Which  multiplies  Cm^, ,  is  an  extremely  small  quantity  when 
compared  with  the  expected  angle  of  attack.  This  is  more  clearly  illus¬ 
trated  by  a  simple  ccryariccr.  of  the  damping  parameter  d.p*  with  the 
static  moment  parcuoeter  v?p“  .  Accordingly,  from  Volume  I,  in  which  the 
satellite  orientation  angles  and  rates  are  expressed  in  terms  of  the 
orbit  parameters. 


-  f  ft r? 
sp  =  ^  ~T  v 
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and,  fron  Eq.  .(172.’. 
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(122) 


023) 


By  inspection  of  Eqs.  (175)  and  (122),  it  :s  seen  that  if  ^  and  P  t  *e 
of  the  same  order  (as  will  be  shown),  then  3*  and  (?/*  differ  in  magni¬ 
tude  by  at  least  the  order  of  c/r p  ;  i.e.,  the  ratio  of  the  satellite 
reference  dimension  to  the  peri gee  radius,  which  for  4  iscnable  satellites 
(witn  dimensions  of  at  most  10?  ft.)  is  of  order  10”5,  *  very  snail  number. 

Although  the  results  of  uie  numerical  analysis  in  Volume  I  adequately 
show  the  insignifi  rsnro  *f  the  aercdyr.am  c  damping  insofar  as  the  satel¬ 
lite's  angular  notion  is  ciKieem^d.  it  is  tne  purpose  of  this  section,  in 
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part,  to  justify  the  original  hypothesis  that  ^  is  not  significantly 
larger  than  V  -  Although  this  investigation  u-d;  first  seen  super¬ 
fluous,  there  ai^  applications  in  which  aerodynamic  damping-m -pitch  may 
conceivably  be  a  useful  input;  e.g.,  the  design  of  a  sensitive  active 
stabilisation  system.  Accordingly,  the  analyses  to  follow  investigate 
the  variations  of  the  dampmg-in-pitch  derivative.  Cm^.  ,  throughout  the 
acute  angle  of  attack  range  for  the  three  basic  body  shapes*  the  cylinder, 
sphere  and  cone.  The  purpose  of  the  angle  of  attack  variation  is  to 
determine  if  there  are  any  large  discrepant  ea  cetween  Cm^  at  angle  of 
attack  and  the  usual  definition  of  Cv*.  ;  i.e.,  in  Eq.  (17o) ,  £►>, ,  is 
defined  as,  for  simplicity,  *  » 


r)c«.  1 


°t=0 


L.2  Danping-in-Pitch  of  Sasic  Sosfr  Shapes 
h.2.1  Introduction 


IXie  to  the  complicated  nature  of  the  integrals,  the  daraping-in- 
pitch  derivative  C**y  is  emulated  using  the  approximate  theory  of 
Ref.  U.  Consequently,  it  m?y  be  shown  that  „,e  evaluation  of  as 

approaches  zero  incurs  very  little  error  since  for  diffuse  reflec¬ 
tion,  the  approximate  theory  yields  a  linear  variation  of  C*v.  with  ^ . 
In  addition,  the  various  F v  ’3  are  calculated  using  optimum  body 
geometry  as  determined  from  the  static  stability  analyses.  Due  to  the 
generally  second-order  importance  of  aerodynamic  damping  as  pertaining  to 
trds  study,  the  optimization  of  body  geometry  for  maximum  aerodynamic 
damping  does  not  seem  justified. 

Diffuse  reflection  1  <r-  cr'-i)  is  assumed,  and  as  before,  the 
reference  area  and  length  are  TTDVd-  and  D  respectively. 


h.2.2  Cylinder 


The  danping-in-pitch  derivative  for  a  cylinder,  without  end  effects 
ana  for  diffuse  reflection,  may  be  evaluated  using  the  approximate  fre*  - 
molecule  flow  theory  (Ref.  u  or  6).  Then,  from  Eq.  (27), 


(178) 
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whe-e  the  cylirder  geometry  is  defined  an  Fig.  31.  Likewise,  the  cmtncu- 
uior  of  ’.he  flat-faced  nose  jorticn  ray  be  f<rna  as 


v'*'WiE 
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bince  mere  is  assumed  to  be  no  interaction  between  the  two  surfaces, 
Sqs.  (1?8)  and  (179)  nay  be  added  to  gave 
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Eq.  (180)  is  presented  for  the  simple  case  of  T w  - T-o  in  Figs.  51,  $2 
and  53  for  cylinder  fineness  ratios  of  1.  2  and  5  respectively.  The 
associated  temperature  correction  term  is  seen  to  be,  from  Eq.  (180),' 


or 


which  is  independent  of  angle  of  attack,  and  is  shown  in  Fig.  5U  fur  a 
wide  range  if  temperature  ratios. 

Unlike  the  case  of  static  pitching  moment  .  ope,  the  non-linear 
variation  of  with  center  of  mass  location  is  clearly  evident.  It 

may  also  be  notea  from  Figs.  5i»  5^  and  53  that  Cn*y.  becomes  more 
negative  (more  damping)  as  angle  of  attack  initially  increases. 


For  purposes  of  evaluating  F y  of  the  cylincer ,  Eq.  (l&O)  for  the 
?ero  angle  of  attack  case  "eccuces  to^ 


x<n 
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The  darping-ir.-pitch  parameter,  ,•  say  now  he  forced  using  Eq.  (132) 

and  Sq.  (12$)  for  the  pitch  ncrccnt  of  inertia.  Hence,, 


(183) 


or  for  the  case  of  an  off-center,  homogeneous  payload,  Eq.  (183)  nay  be 
written  in  terras  of  payload  volume  ratio  anc  weight  using  Eq.  (127)  and 


the  assumption  that  m fv*  j  !.€.. 

r  ~  /  /  T-  S  fi-)V  J f-  5./3C  f 
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*It  ic  real.-  od  that  the  "shadow”  area  on  a  cylinder  may  have  a  signifi¬ 
cant  contribution  to  Ci*j  for  c<  near  zero.  For  the  practical  case  of 
low  *(D  and  hign  ,  however,  this  contribution  does  not  alter 
Sq.  (182)  appreciably. 
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For  a  sample  calculation,  let  <'’r/'*r  equal  1/2  and  3/a,  and  let  the 
respective  fineness  ratios  be  1.80  ana  1.08,  as  ta»cen  from  fable  1.  Then 
Kq.  (18k)  becoMf,  ror  (l/S„)  [t^.  S  .  i45  , 
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Plots  of  3qs.  (3  85)  are  pi  esex«ted  in  Figs.  61  and  62  along  with  similar 
vax iato.o.is  of  \  y  *or  ^r,e  sphere  and  cone.  Discussion  of  these  rigures 
is  reserved  for  Section  U-2-5- 


h.2,3  Sphere 


The  approximate  damping-in-pitch  derivative  of  a  sphere  nay  be 
found  by  integrating  Sq.  (27)  ever  t<ie  front  (or  windward)  hemisphere. 
For  diffuse  reflection,  the  result  is  readily  found  to  be 


Eq.  (1S6)  is  shown  in  Fig.  55  for  T„ 
correction  tern,  wnich  is 


t=t)  +«‘*]{,-;Tf}  •  <i56> 


-=T„ 


The  normalized  temperature 


A- 


(187) 


is  plotted  in  Fig.  $6.  Unlike  the  cylinder,  the  sphere  damping  is  seen 
to  decrease  with  increasing  angle  of  attack. 

Evaluating  Eq.  (186)  at  zero  angle  of  attack  loads  to  the  usual 
definition  of  , 


7.  ijQL  ifC  /. .  z’Jaj-'  ,V !  -z^  . 
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For  -ho  case  ci‘  a  "loadea"  3phexe  with  a  horiogeneois  p^ylv-a'..  substi¬ 
tution  of  Sq.  (138)  and  ?/;.  (133),  into  2u,  (176;  results  in  an  expression 
lor  the  sphere  ;  i.e., 


,  r  «*  * ,  _  Xfu/i  i  / 1  *>  v 

r,  =  irD^B-z^lTg.  )  V.^‘-2tJ  +  jU~  2  u  ) 

9»'t.  f  >  -  ii-£  +  k-e’-  -  If’ 
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The  Centex  of' mass  nzy  oe  written  in  terms  of  depth  ra*io  t  ,  and  the 
aia.-r.eter  nay  be  expressed  as  a  function  of  payl*~j  vsight  foz  a  specified 
reload  density.  Then  by  using  Fqs.  (13u )  and  (lj5)  and  r./v  s 

fed «/«j>  £<i-  (iSs,)  bf'coT-°s 
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where  £  is  a  function  of  tne  payload  volume  ratio  f roM  the  first  of 
Eqs.  (135)  or  Fig.  33-  For  the  sanple  volume  ratios  of  1/2  and  )/U,  tne 
sphere  is  approximately 
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where  C1/^)  nV/T«©  is  lanen  to  be  .065.  Tne  foregoing  expressions 
for  are  plotted  m  Figs.  61  and  62.  A  discussion  of  tnese  results 

is  2i  Section  h.2.5- 


L.2.1  Cone 

Following  the  example  of  cone  rtatic  pitching  ^  '■-nt,  tne  Hawping - 
in-pitch  derivative  C*,y  may  be  estimated  using  the  same  approximate 
theory.  Then,  for  njffdsa  reflection,  the  it.tegrovicn  of  (27 )  yields 
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where  4  H  W/+<*«.<  .  as  before,  ana 

mW 

=  10  [4  1 1  c] 
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Ccr.requently,  for  (1??)  «r.plilles  to 
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and  for  o<  identically  equal  to  zero, 
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or  re  awwSMST?11”*  \’h8r«rt  “  *“‘e  ^.isate  theory,  a  plot 

A  S;  S2  0°  f  *?•  e-  at'“k  r“?e  «M  shop  a  sharp  !U*i: 

r  ,  ?  a  y"  sus  «  at  the  points  where  c<  o  £  ,  vhe-r 
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ho  fa.rly  accurate  for  the  calculation  of  cone  lift.  dr“r  and 
of  'cj1  &JTT’  “^t5a?t  is  "**  "troin  to  e™0th  trie  variation 
«nti^-™hi  ^  :slier -nf  thc  •trteta  of  the  »nadou.  This  refine- 

-o’ecu’a^-Me'^r  t?  y  U”a  ■3f'J'’'ir'K  *"  avera;'<‘  nor°al  coaponeot  of 

“Peed  '“u°  ov,,r  the  s.'.sdovj  and  than  substituting  th<n  cunnt.<‘- 

[?frror  £?£$?  and°f inall  v^i *  “*T”“  r0r  C~V  (derivative "of he 
>  ^  '*  finally  integrating  over  the  shadow  portion  Thp 
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written  as 
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Without  elaboration,  the  result  is 
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In  Eq.  (197),  is  expressed  by  Eq.  (22),  oy  Eq.  (hi)  and 

p(<)  by  Eq.  (h2).  Also,  the  shadow  angle  is  equal  to  cos*'(-<)  = 

TT  PC^)  - 

The  final  cone  G*y  equations  -  Eq.  (19U)  for  oC  ^  S  and 
Eq.  (196)  for  o'  i  9>  -  -re  plotted  in  57,  58  and  59  for  cone 

angles  cf  15°,  30°  and  15°  respectively.  As  before,  Tw  is  set  equal 
to  T»  for  simplicity.  Fig.  60  presents  the  corresponding  variation 
cf  AC*n^  (normalized)  due  tc  Tw  different  iron  T*,  .  The  equations 
xor  these  curves  are 
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It  say  be  seen  in  Figs.  57  through  60  that  the  breaks  in  the  curves  at 
oc  =  £  are  still  evident,  especially  m  Fig-  60,  and  for  $o=°o  lr. 
Figs.  57  through  59.  The  use  of  the  shadow  corrt.wticn,  however,  defin¬ 
itely  improves  the  smoothness  of  the  Cmy  curves  at  the  lower  speed 
ratios.  * 
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New  that  has  oeen  formulated,  the  damping;- m-pitch  parameter 
*V  may  he  written  using  Sq.  (195)  and  Eq.  (139)*  Accordingly,  Ec .  (176) 
becomes 
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Using  the  optimum  cone  angle,  S  =■  o>f‘  fzT  .  from  Section  3.U.2.3,  and 
the  relation,  x«v*/ JL  -  (3/4-)(nr/v-)  v3  ,  Eq.  (199)  may  be  expressed  in 

♦eras  of  the  payload  volume  rr tic .  Hence, 
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•-here  f ff5  .  For  the  Fresc-nt  examples.'  the  substitution 

of  •v'/or-  =  i/?  and  3A  into  Eq.  (200)  yields 
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Piets  cf  Eqs.  (2“Ji)  are  ahown  in  rigs.  61  and  62  along  with  the  corres¬ 
ponding  cylinder  and  sphere  expressions. 
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h .2 .5  Comparison  of  Basic  Body  Snapes 

As  noted  previously,  Kigs.  61  and  62  present  plots  of  the  dynamic 
s  tability  parameter  Py  versus  satellite  payload  might  for  the  special 
cases  ’-here  the  “loaded**  'onfiguraticr.u  hs**c  payload  volume  ratios,  of  1/2 
ar.tf  j/u.  The  p^yloau  is  assumed  homogeneous  with  a  mass  density  of 
l/?  slug/ft'.  The  figures  clearly  indicate,  as  expected,  that  the  cone 
r»as  more  aerodynamic  damping  thar.  cithei  wie  cylinder  or  sphere;  at  least 
for  the  optimum  geometry  (for  maximum  f  )  assumed  for  the ;e  examples. 

Ely  comparing  rigs.  61  and  62  with  rigs.  36  and  37,  it  is  aiso  noticed,  as 
hypothesised  previously,  that  Pv  13  of  the  same  crcer  of  magnitude  as 
the  static  parameter  P  ,  especially  at  the  lower  weights.  Consequently, 
it  ray  b**  concluded  that,  since  Py  is  mot  significantly  larger  than  P 
(say  by  a  factor  of  lO1),  the  assertion  in  Volume  I  of  essentially  negli¬ 
gible  aerodynamic  darling  may  be  regarded  as  fully  substantiated. 


Cylinder  Datipln8-ii-Pitch  Oerivati, 
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CONCLUSIONS 


Recent  dev-lojr.ents  in  the  application  of  free- molecule  theory  (Ref.  6) 
have  permitted  a  rather  detailed  aerodynamic  analysis  of  near-Earth  satel¬ 
lites.  This  theory,  although  restricted  to  convex  body  shapes,  has  been 
utilised  tc  yield  rapid  and  nearly  exact  aerodynamic  force  and  moment  calcu¬ 
lations  for  most  satellite  configurations  and  orbit  conditions.  Moreover, 
tnere  are  no  restrictions  as  to  the  variance  of  either  the  tangential  or 
normal  moff<entum  exchange  coefficients,  5*  and  C*'  .  As  shown,  the  surface 
temperature,  Tw  ,  may  also  be  treated  as  a  variable  over  the  body.  In  view 
vt  the  general  nature  of  this  approximate  theory,  numerous  calculations 
(wr.ich  riav  be  too  detailed  for  some  readers)  and  charts  have  oeen  presented 
33  examples  of  the  applicability  of  the  theory.  Although  the  foregoing 
Study  was  intended  to  "stand  alone",  tr.e  primary  purpose  of  tins  report  was 
to  augment  the  conclusions  of  Volume  I  of  this  series  of  reports.  That  is, 
aerodynamic  stabilization  and  orientation  of  near-Earth  satellites  is 
feasible  as  long  as  proper  (but  not  ’inreas enable)  aerodynamic  design  and 
internal  mass  distribution  conditions  are  met. 

A  supplemental  but  significant  part  of  the  stud?  was  concerned  with 
estimating  drag  coefficients  of  simple  configurations,  ine  calculations 
showed  that  if  the  drag  coefficients  are  based  on  the  projected  frontal 
area  (quite  often  a  random  variation)  of  the  satellite,  a  unified  drag 
coefficient  variation  could  be  approximated  far  v-ost,  lo«  fineness  ratic 
(order  of  one)  configurations.  For  long  cylinders  or  very  sharp  cones, 
however,  the  quasi-unucrs  drag  variation  with  molecular  speed  ratio  could 
be  in  considerable  error. 

Since  the  primary  purpose  of  this  study  was  to  investigate  the  effects 
of  aerodynamic  torques,  subsequent  portions  of  this  volume  were  concerned 
with  methods  for  estimating  aerodynamic  static  and  dynamic  pitching  moment 
coefficients  as  well  as  the  incorporation  of  these  coefficients  in  the 
design  of  practical  satellite  configurations.  Although  the  calculation  of 
aerodynamic  coefficients  in  fre  '•-molecule  flow  are  inherently  different  from 
that  in  continuum  flew,  the  conventional  "airplane"  definitions  of  C**  . 

Cm*  Cinv  were  conveniently  utilized  throughout  this  volume.  Expres¬ 
sion?  for  the  static  pitching  moment  coefficient.  Cm,  ,  were  presented  for 
the  three  basic  cylindrical,  spherical  and  conical  body  shapes.  In  addition, 
charts  of  Cm,  versus  angle  of  attack  were  presented  for  these  shax>es,  whit 
may  be  useful  in  other  free-molecuie  design  studies  not  related  to  this 
report.  These  charts  were  also  uccd  tc  substantiate  the  linear  approxi¬ 
mation  for  **r*-»ll  angle  of  attack  ranges j  a  basic  and  com*.  assumption. 
Consequently,  expressions  for  the  initial  pitching  moment  slope.  Cm*  ,  were 
formulated  for  the  basic  body  shapes,  and  suosequently  for  the  more  complex 
shapes  of  the  power  law  body,  the  truncated  cone,  the  hemisphere -cylinder, 
and  the  eor.e-e/linder.  For  simplicity  purposes,  the  usual  assumption  of 
diffuse  molecular  reflection  was  utilized  for  these  Cm*  calculations. 
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Since  the  free-r'o'lscule  theory  used  in  this  study  was  not  restricted 
to  specific  values  of  the  i  effect  ion  coefficients,  <r  and  <r'  ,  and  the 
surface  terete ratures,  Tv  ;  examples  were  presented  whicn  illustrated  the 
effects  of  oo  tli  uniform  and  non-uni  form  vaiiations  of  (T  ,  and  Tw  . 

It  was  shown  that  the  surface  temperature  level  as  well  as  longitudinal 
and  axial  variations  of  temperature  have  little  effect  on  static  stability 
or  orientation,  at  least  for  the  simple  configurations  considered.  On  the 
other  hand,  the  magnitude  of  the  reflection  coefficients  was  shown  to 
affect  the  pitching  moment  slope  by  varying  amounts,  depending  on  the 
configuration  and  the  relative  dominance  between  pressure  and  shear,  for 
nearly  diffuse  molecule-surface  interactions  (  <T  and  0*‘  between  0.9 
end  1.0),  however,  the  example  pitching  moment  slopes  were  shown  to  be 
only  slightly  different  from  completely  diffuse  values. 

Based  upon  the  satellite  motion  analysis  in  Volume  2,  txie  aerodynamic 
stability  parameter  V  was  investigated  for  the  three  basic  body  shapes 
as  veil  as  the  more  complex  configurations.  Since  P  was  defined  as 
being  the  pitching  moment  times  its  reference  dimensions,  divided 

by  the  pitching  moment  of  inertia,  the  maximization  of  P  established 
certain  optimum  body  geometry  for  each  set  of  constraints.  The  satellite 
models  selected  for  the  investigations  had  homogeneous  payloads  concen¬ 
trated  in  the  forward  portion  with  rigid,  but  massless,  shells  cr  skirts 
comprising  the  aft  portions.  The  skirts,  which  could  logically  be  fairings 
between  the  payload  and  booster,  augmented  tne  aerodynamic  stability  with¬ 
out  affecting  the  pitching  moment  of  inertia,  'ccordingly,  the  optimiza¬ 
tion  process  tended  to  maximize  the  aerodynamic  contribution  while 
minimizing  the  mass  contribution.  For  a  specified  payload  volume  to  total 
satellite  volume  ratio  and  a  specified  payload  mass  density,  the  maximiza¬ 
tion  of  P  limited  in  optimum  configurations  vhith  had  rather  low  (and 
practical)  lergth  to  diameter  ratios.  On  this  basis,  it  was  shown  that 
the  pointed  core,  the  power  law  body,  and  the  truncated  cone  were  generally 
the  most  stable  for  a  given  payload  weight.  For  other  constraints,  however, 
such  as  a  specified  maximum  diameter,  the  cylinder  and  the  cone-cylinder 
showed  precise  ju»  being  useful  and  practical  aerodynamic  sbipoc.  Tho 
sphere  with  its  impractical  s»^?il  was,  as  expected,  the  poorest  configura¬ 
tion.  The  optimum  sphere-capped  cylinder,  however,  was  shown  to  have 
nearly  t'ne  some  static  stability  as  the  optimum  flat-faced  cylinder. 

Since  the  approximate  free-molecule  theory  used  in  this  study  permitted 
rapid  calculations  of  the  aerodynamic  damping-in-pitch  derivative,  and  si  -e 
pitch  damning  was  necessarily  considered  in  the  satellite  motion  analysis 
(volume  I),  the  aerodynamic  investigation  vas  conclude  with  a  brief  study 
of  the  daming-in-pitch  derivative,  Cmy  >  and  its  cou.  Verpart,  Py  ,  <*or 
the  besic  cylinder,  sphere  and  cone  configurations.  As  in  the  case  of 
static  pitching  moment,  C>*.  ,  design  charts  of  were  presented  for 

the  basic  body  shapes.  Since  it  was  argued  in  Volume  I  that  aerodynamic 
damping  is  negligible  if  Py  is  the  same  order  of  magnitude  as  P  , 
e»l ml aMo«s  of  f—  wye  made  for  the  cone,  sphere  and  cylinder 
-isir.tf  oplimcft.  static  stability  geometry  (for  maximum  V  ).  The  insignifi¬ 
cance  of  aerodynamic  damping,  as  pertaining  to  this  study,  was  verified. 


Although  this  study  has  provided  a  foundation  for  sateliitu  aero¬ 
dynamic  desigr-,  several  important  and  unansw»rec‘  questions  remain.  For 
examble,  the  lack  o’'  a  general,  three-dimensional  tneory  for  predicting 
the  effects  of  molecular  reflections  between  surfaces  has  prevented  t*c 
aerodynamic  analysis  of  concave  bodies;  3.g.,  configurations  with  flared 
shirts  cr  inflatable,  spherical  stabilisers.  Trac  lack  of  a  usable 
surface-surface  interaction  theory  has  also  excluded  investigations  of 
ccreplex,  "winged"  satellite  configurations,  especially  from  the  viewpoint 
of  orientation  and  orbit  control.  In  addition,  aerodynamic  analyses  of 
satellites  and  sacelloids  in  the  near  free -molecule  or  transitional 
regime  <ue  essentially  non-existent.  Consequently,  it  is  clear  that  the 
practical  applications  of  rarefied  gasd/rumics  have  been  only  briefly 
exploitea.  Sven  with  the  extensile  analysis  contained  Kerein,  the  free¬ 
rider  vie  aerodynamic  investigation  of  simple,  convex  body  shapes  is  by 
no  means  complete.  For  instance,  the  simplifying  assumptions  of  homo¬ 
geneous  payloads  and  massless  skirts  leave  much  to  be  desired.  Also,  no 
attempt  was  made  to  analytically  determine  an  optimum  configuration,  by 
the  methods  of  variational  calculus  techniques,  which  would  provide  an 
absolute  maximum  for  the  static  stability  parameter  f  . 

Insofar  as  passive  aerodynamic  stabilization  and  orientation  of 
near -Earth  satellites  has  been  foind  to  be  feasible,  it  is  recommended 
that  concentrated  effort  be  given  to  the  advancement  of  applied  rarefied 
gasqynamics  in  order  to  lay  the  foundations  for  active  aerodynamic 
satellite  control.  It  is  realized,  of  course,  that  a  conceptual  under¬ 
standing  of  basic  gas  kinetics  is  a  prerequisite  to  the  development  of 
a  useful  and  accurate  method  of  predicting  lv*  density,  aerodynamic 
characteristics.  This  also  requires  increased  study  of  the  molecule- 
surface  interaction  phenomenon,  especially  at  satellite  energies. 
Accordingly,  numerous  problems  remain  which  both  research  and  engineer 
rust  solve  if  aerodynamic  control  and  stability  is  to  be  incorporated  in 
sophisticated  satellite  systems. 
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7HF  EFFECT  Or  SOLAR  RADIATION  PRESSURE  ON  SATELLITE  TORQUE 


Due  to  tr*e  (nearly)  exponential  decrease  of  atmosphere  density  with 
altitude,  indications  are  that  the  magnitude  of  satellite  torques  aue  to 
solar  radiation  pressure  roay  approach  the  magnitude  of  aerodynamic  ana 
gravitational  torques  at  altitudes  i.eai  hw  miles.  In  oi u«r  to  estimate 
these  solar  torques  more  precisely,  it  is  r.ecessarv  to  develop  expressions 
for  photon  pressure  and  shearing  forces  cn  a  body  surface  of  genet al  shape. 
The  qualitative  analysis  tu  follow  assumes  that  the  free -stream  photon 
Bgas  a*?  analogous  to  a  unidirectional  particle  field,  ar.j  that  the 
photons  have  a  mass  density  and  travel  with  the  speea  of  light.  Accord¬ 
ingly,  this  analogy  permits  utilization  (out  slight  modification;  of  trie 
fr ee-nolecule  aerodynamic  expressions  in  Section  1. 

The  pressure  and  shearing  stresses  due  to  an  oncoming  photon  stream 
Incident  upon  a  unit  surface  area  nay  be  respectively  written  as 


p. 


'—CD 


(A-l) 


where  Ps^  repz esents  the  free-stream  radiation  pressure  (one-half  of 
the  maximum  measurable  pressure)  and  where  c*j  and  Ct  are  the  normal 
and  tangential  components  of  the  free-stream  photon  velocity  C*»  .  If 
all  the  incident  photons  arc  reflected  from  the  surface  specularly,  then 
normal  and  tangential  forces  per  unit  area  due  te  the  rebounding  photons 
are 


Pr=Pi  -  IV.  g. 

'V  •=  p  CmCj 

CZ 


Ho  absorption ,  | 

specular  reflection] 


(A-2) 


and  her.ce,  the  total  force  components  per  unit  a.  ia  become 


Vi' 


P=  pt-pr  = 

V  s  'Y:  s  o 


which  arc  analogous  to  the  Newtonian  equations  used  in  aercdyrim.cs  if 
P<  -x,  is  replaced  by  the  free-stream  dynamic  pressure.  Eqs.  (\-3) , 
iiv>«)vvai  |  do  ijv  v  repieser.t  the  photon  reflection  process  existing 

or.  s.ost  satellite  materials.  That  is,  (l)  sorce  of  the  photons  are  not 
reflected  but  may  be  absorbed  and  transmitted  through  the  material,  and 
(2)  the  reflected  photons  iaay  not  rct^^r.d  in  a  specular  fashion.  Ar. 
extrere  of  the  latter  is  the  case  of  completely  diffuse  reflection  in 
which  there  is  no  preferred  direction  for  the  rebounding  pnotons.  Then 
if  the  actual  reflection  process  is  assumed  to  vary  linearly  between 
specular  and  diffuse,  fiqs.  (A-2)  may  be  more  correctly  written  as 


Pr  =  a[(!-A)0,  +A  Pj] 
Tr  =A[(l_A;'f'i  4  ATjl 


(A-ll) 


where  the  subscript  d  refers  to  completely  diffuse  reflection;  the 
reflectivity  A  represents  the  percentage  of  incident  photons  reflected; 
and  A  is  a  erasure  of  the  diffusivitr:  i.e.,  i  is  for  completely 
diffuse  reflection  while  A  =  O  is  for  completely  specular  reflection. 
Since  there  is  no  preferred  direction  for  diffuse  reflection.  *Yj  2  O 
Furthermore, 


Pf  3  I'-i  '  ' 

Tr  =  A.  (I  -  Ay  7- 


(A-5) 


and  the  total  pressure  and  shearing  stress  may  be  written  as 


P  =  P,  +  Pr  ■={  l  + 

f  1 

r  -  'V;  - yr  =  ! !  -  a  4  ax J  . 

L  Cco 

For  ccapl«*tely  specular  reflection,  a  -  O  ,  and  Eqs.  (A-o)  oeconc 


P= 

Coo 

r  =  C-^|T 

Coo 

Likewise,  lor  diffuse  reflection,  X  =  i  ,  and 


D  =  (  |  +  4*)ps 

‘  '  S  /KS"C>  (A-8 

r  =  p<  . 

1  *  f  T- 

For  the  examples  to  follo<rj  i.e.,  the  effect  of  solar  radiation 
pressure  on  satellite  ♦/'«*^'iec,  the  specular  equations.  Eas.  (A-7)  will  be 
assuneu  for  siraplicity . 

>.l ih  the  foregoing  equations,  the  solar  torque  on  a  body  of  revolu¬ 
tion  xay  oe  calculated  using  the  last  of  Sqs.  (20),  which  is  rewritten 
nere  as. 


+  rf  — -Y  i  +  r'1  <•(*  ■  Xu^  r  <)-..<* 


(A-?) 


is  th- 


ViS 


wrier  a  the  integration  is  ove»  the  sunny  surface  only  3nd  wnere  <^5 
•solar  'angle  ol  atta*°-'T.  Aiso,  from  Sq.  (22), 


Con 
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For  the  present  study,  it  is  of  more  interest  to  calculate  the  initial 
torque  slope  (i.e.,  a  measure  oi  the  staoili*y)  ratncr  tnan  tnc  total 
torque.  Then  by  comparing  tne  solar  torque  slope  with  the  aerodynamic  and 
gravitational  torque  slopeo,  the  “homing*1  cnarac  tori  sties  of  a  satellite 
rsay  oe  evaluated.  S^nce  the  cylinder,  sphere,  and  cone  have  been  exten¬ 
sively  treated  m  the  aerodynamic  ar^lysis,  it  is  eonven?  ent.  to  calculate 
the  soiar  torque  slopes  for  th»-3e  bodies.  Thus,  the  cylinder,  sphere  and 
cone  evaluations  of  the  derivative  of  Eq.  (A-/)  with  respect,  to  o'*  in 
the  limit  — y  o  are 
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wncre  cylinder  nose  effects  are  included  and  vhe~e  the  specula  phctwi 
reflection  expressions,  Eas.  (A-7),  are  used. 

The  corresponding  aerodynamic  torque  slopes  are,  from  Eqs.  (6u), 
(6?)  and  l?3)# 


1ii6 


i 


'  , „ 


*  '  C*M  u  * 


iin 

[>  i  1.0  J  u* 


,  !  f  f 

*'ys- 


i  XcvrJ^ 


/  }u\ 

i~5T) 

V  ■'“'yiPntaf. 


4-11) 


/iw 

\  n-  ) 

'  •  *  >  Out 


+ 


where  ^TO  is  the  free -stream  dynamo  pressure. 

Finally,  the  gravitational  torque  slopes  nay  be  found  from  Appendix  A 
of  Volume  I.  For  a  body  cf  revolution,  the  initial  torque  slope  s 


wnere 
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/f  »  product  of  Earth* s  nass^and  universal  agitational 
constant  =  l.uuo  x  it^/eec*- 

r  ■  distance  from  Sarth*s  center  to  satellite  mass  center 

€  *  inclination  of  the  satellite  X  axis  to  the  local  horizon. 


7hs  moments  of  ine:  tia  of  the  cylinder,  sphere  anc  cone  havt  teen  deter- 
rained  in  Section  3.U.2.  Consequently,  Eq.  (A-12)  becomes 
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~r.cre  e.  is  the  sphere  payload  fineness  ratio  aid  a  function  of  Xunsl  £5  . 

As  an  example ,  Eqs.  (A-10) (A-llj  and  {A -13)  are  compared 
Figs.  A-l.  A-2  and  A- 3  for  the  special  case  of 

payload  volume,  'sr  -  10  ft^ 
payload  volume  ratio,  />rl'v r  ■  l/2 
payload  density,  m/'v-  •  j/2  . 

The  optimum  aerodynamic  geometry  from  Section  }M. 2  is  selected  for  these 
conditions.  Also,  a  circular  orbit  is  assumed  so  that  o%«  aay  be  written 
in  terms  or  altitude  (Kef.  7).  The  reflectivity  /v  is  taken  to  be  0.5, 
and  the  free-stream  solar  pressure  Pic*  is  assumed  constant  nr«d  *oua1 
to  10””*  Ibs/ft*  (an  average  of  Refs.  Hi,  15>  and  16)..  As  before,  Tw  is 
assured  to  be  550°?.  an-4  hence,  Fir.  2  is  applicable. 

Tne  curves  in  Figs.  A-l,  A-2  ai.d  A-3  clearly  show,  at  least  for  these 
sample  configurations,  that  solar  torques  are  small  vnen  compared  with 
aerodynamic  torques  for  altitudes  below  2  million  feet  (say).  Moreover, 
if  the  aerodynamic  and  gravitational  torque  slopes  are  combined  (a  'sacure 
of  <3f*  ).  J  ien  the  importance  oi  the  gravity  gradient  at  the  higher 
altitudes  is  clearly  illustrated.  That  is,  the  cohere  and  core  (Figs.  A-2 
ana  *.  })  have  stable  gravitational  torque  slopes  at  ^  -  o  /t  and  hence, 
augaent  the  aerodynamic  stability.  Accordingly,  the  dominance  of  solar 
torques  is  delayed  until  higher  altiu*dcc  are  reached.  The  opposite  effect 
is  evident  in  the  case  of  the  cylinder  (Fig.  A-l)  which  has  an  unstable 
gravitational  orq„c  slope  (for  this  example);  i.a.,  the  solar  torque 
becur.es  dominant  at  a  lower  altitude  due  to  the  reduction  of  the  combined 
aerodynamic  and  >ravitatioral  stability. 
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